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1.0  INTRODUCTION 

The  transmissions  of  extremely  low  frequency  (ELF)  waves  at  nighttime  are 
often  subject  to  amplitude  attenuations  of  several  dB  (Davis,  1974,  1976; 

Davis  and  Meyers,  1975;  Larsen,  1974).  A  successful  correlation 
between  anomalous  signal  strengths  and  precipitating  particles  was  established 
for  the  first  time  in  1976  (Imhof  et  al,  1976;  Reagan  et  al,  1978  ).  This 
preliminary  finding,  based  primarily  on  the  measured  transmissions  between 
the  Wisconsin  Test  Facility  (WTF)  and  the  mid-latitude  receiving  site  in 
Connecticut  and  direct  satellite  measurements  of  precipitating  particles, 
has  formed  a  basis  for  further  more  quantitative  studies  of  the  cause-and- 
effect  relations!. ip. 

Subsequent  studies  of  the  effects  of  energetic  particle  precipitation 
on  ELF  transmission  signal  strength  have  been  made  at  the  Lockheed  Palo  Alto 
Research  Laboratory  (LPARL)  based  partly  on  data  acquired  with  scientific 
payloads  on  the  low  altitude  polar  orbiting  satellites  1971-089A  and  1972-076B 
developed  by  the  Space  Sciences  Laboratory  of  LPARL  for  the  Office  of  Naval 
Research,  the  Defense  Nuclear  Agency,  and  the  Defense  Advanced  Research 
Projects  Agency.  In  these  investigations,  supported  by  the  Office  of  Naval 
Research,  ELF  propagation  data  were  obtained  from  Dr.  John  Davis  of  the  Naval 
Research  Laboratory,  Washington,  D.C.,  and  P.  R.  Bannister  of  Naval  Underwater 
Systems  Center,  New  London,  Connecticut.  This  large  data  set  includes 
measurements  taken  during  coordinated  exercises  involving  ELF  transmission 
between  the  U.S.  Navy  Wisconsin  Test  Facility  and  receiving  stations  in 
Connecticut,  Maryland,  Greenland,  Norway,  and  Italy  performed  in  March  -  April 
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1975  and  March  -  April  1976.  These  studies  have  provided  further  verification 
of  the  importance  of  particle  precipitation  and  have  also  demonstrated  the 
influence  of  sporadic  E-layers  on  ELF  signal  strength.  Since  the  details  of 
the  horizontal  and  vertical  distributions  of  the  enhanced  ionization  are  clearly 
very  important,  a  large  body  of  coordinated  data  is  required  to  provide  quan¬ 
titative  interpretations  of  the  results  and  to  assess  the  impact  of  important 
chemistry  parameters  on  the  transmitted  ELF  signal  strength.  Clearly,  good 
chemistry  models  are  also  needed  for  a  proper  evaluation  of  the  coordinated 
data.  In  this  regard  the  satellite/Chatanika  backscatter  radar  data  acquired 
during  the  August  1972  solar  particle  event  have  been  particularly  useful. 

In  1976  a  qualitative  correlation  was  established  between  anomalous  ELF 
signal  levels  received  at  Connecticut  and  the  precipitation  of  significant 
fluxes  of  electrons  from  the  radiation  belts  as  measured  with  the  LPARL  pay- 
loads  on  the  low  altitude  polar-orbiting  satellites  1971-089A  and  1972-076B. 
However,  a  detailed  quantitative  correlation  between  the  sign  (signal  enhance¬ 
ment  versus  degradation)  and  severity  of  the  anomalies  and  particle  character¬ 
istics  such  as  intensities  and  spectra  was  not  firmly  established  because  of 
the  limited  data  base  used.  Stimulated  by  these  findings,  special  coordinated 
exercises  were  conducted  in  March  -  April  1976  involving  satellite  measurements 
of  the  precipitating  particles  and  ELF  transmissions  between  the  U.  S.  Navy 
Wisconsin  Test  Facility  and  receiving  stations  in  Maryland,  Greenland,  Norway, 
and  Italy.  Coordinated  satellite/ELF  transmission  measurements  were  also  made 
in  March  -  April  1975  and  July  1975. 

From  coordinated  satellite  and  ELF  attenuation  measurements,  it  has  been 
found  that  direct  particle  precipitation  into  the  atmosphere  can  cause  ELF 
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transmission  anomalies.  In  these  anomalies  the  signal  strengths  may  be  either 
attenuated  or  enhanced  depending  upon  the  geometry  and  details  of  the  ion  and 
electron  density  profiles  resulting  from  the  particle  precipitation.  Detailed 
discussions  of  the  investigations  performed  from  the  coordinated  data  sets 
acquired  in  1975  and  1976  are  provided  in  the  previous  annual  reports  (Imhof, 
et  al.,  1976,  1977,  1978)  and  by  Reagan,  et  al.  (1978). 

The  purpose  of  the  present  report  is  to  present  further  work  on  the 
influence  of  precipitating  energetic  particles  on  ELF  propagation.  The  effort 
principally  involves  the  effects  of  solar  particle  events  in  enhancing  the 
conductivities  of  the  earth's  ionosphere  wave  guide.  Calculations  have  been 
made  of  the  sensitivity  of  the  ELF  signal  strength  to  changes  in  the  conduc¬ 
tivities  over  the  transmitter  that  may  occur  during  SPE  conditions.  In  order 
to  be  able  to  calculate  more  accurately  the  changes  in  ELF  signal  strength 
associated  with  solar  particle  events  and  other  disturbed  conditions  at  twilight, 
new  electron  density  data  measured  in  the  4  August  1972  event  have  been  used  to 
provide  increased  knowledge  of  twilight  behavior  which  will  lead  to  improvements 
in  ion  chemistry  modelling.  The  effects  of  the  charge  deposition  in  the  polar 
caps  on  the  earth's  electric  fields  and  current  systems  have  been  modeled. 

More  realistic  ambient  electron  and  ion  density  daytime  profiles  have  been 
proposed  to  provide  accurate  baseline  signal  strength  calculations. 

The  changes  in  conductivity  at  times  of  SPE  events  discussed  above  have 
all  been  based  on  calculations  of  the  effects  of  the  measured  energetic 
particle  inputs.  It  has  been  recognized  that  these  calculated  effects  should 
be  confirmed  experimentally  in  a  VLF-ELF  field  program.  Accordingly,  LPARL  has 
participated  in  appropriate  planning  discussions  and  further  investigation 
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has  been  made  of  an  acceptable  criterion  for  deciding  if  a  SPE  event  is  suitable 
for  an  ELF  propagation  field  experiment  based  on  satellite  measurements  of  pre¬ 
cipitating  protons. 
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2.0  PROPAGATION  OF  ELF  WAVES  -  THEORETICAL  APPROACHES  TO  THE  CALCULATION 
OF  FIELD  STRENGTHS 

The  propagation  of  ELF  waves  is  subject  to  anomalies,  some  of  which  are 
believed  to  be  of  ionospheric  origin,  but  the  causes  and  effects  are  not  well 
understood  at  present.  These  anomalies  result  both  from  signal  absorption  and 
from  changes  in  excitation  of  the  waves  into  the  propagation  medium,  a  wave¬ 
guide  formed  by  the  earth  and  the  lower  ionosphere.  The  conditions  of  the 
ionosphere  are  subject  to  major  day/night  differences  and  indeed  the  ELF  signal 
strengths  do  display  diurnal  variations.  Calculations  of  ELF  transmissions 
have  shown  that  the  signal  strengths  may  depend  in  a  complex  manner  on  the 
ionization  profiles  in  the  waveguide.  The  conductivity  of  the  earth  depends 
on  whether  the  propagation  is  over  rock,  sea  water,  ice  cap,  etc.  The  conduc¬ 
tivity  of  the  atmosphere  depends  on  the  deposition  of  ionizing  particles  in  the 
atmosphere  and  the  absorption  of  ultra  violet  radiation  from  the  sun.  Below 
40  km  altitude,  the  conductivity  of  the  atmosphere  is  produced  by  the  deposition 
of  energy  from  galactic  cosmic  rays  during  quiet  times,  and  by  the  deposition 
of  energy  from  solar  cosmic  ray  protons  during  solar  particle  events.  Unfor¬ 
tunately,  there  has  been  little  experimental  verification  of  the  methods  used 
to  determine  the  conductivities.  At  altitudes  above  ^  50  km  the  electron  con¬ 
ductivities  are  affected  by  electron  precipitation  but  they  can  be  predicted 
from  the  ion  production  rates  with  reasonable  accuracy.  In  strong  events  the 
Bremsstrahlung  x-rays  produced  by  precipitating  electrons  may  cause  significant 
ionization  at  lower  altitudes. 

In  this  chapter  we  will  briefly  describe  some  of  the  approaches  that 
have  been  used  to  calculate  ELF  field  strengths. 


2-1 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
iockmccd  Minim  a 


s  p  a  c  i 


COMPANY 


N  C 


LMSC  D767238 


2 . 1  The  NOSC  Model  Computer  Program 

A  waveguide  model  computer  program  extensively  used  over  the  last 
years  and  which  is  also  being  used  by  the  present  authors,  was  developed 
at  Naval  Ocean  Systems  Center  (NOSC)  and  is  well  documented  in  the  literature 
(see  e.g.  Sheddy,  1968;  Pappert,  1970;  Pappert  and  Moler,  1974). 

The  program  is  based  upon  the  formalism  developed  by  Budden  (1961) 
and  Budden  and  Daniell  (1965)  for  electromagnetic  wave  propagation  in  the 
terrestrial  waveguide.  Special  emphasis  has  been  given  to  its  applicability 
in  the  extremely  low  frequency  (ELF)  range. 

The  program  allows  for  a  terrestrial  waveguide  with  the  following 
characteristics : 

-  Arbitrary  electron  and  ion  density  distributions  with  height 

-  Variable  electron  and  ion  collision  frequency  height  profiles 

-  The  lower  boundary  of  the  waveguide  is  assumed  to  be  a  smooth 
homogeneous  earth  with  arbitrary  conductivity  and  permittivity 

-  Allowance  is  made  for  the  anisotropy  of  the  ionosphere;  arbitrary 
values  for  the  strength  and  direction  of  the  earth's  magnetic  field 
may  be  applied 

-  Horizontal  inhomogenet ies  along  the  propagation  path  may  be  taken 
into  account  by  segmenting  the  path  into  regions  having  different 
values  for  the  waveguide  characteristics  listed  above. 
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The  computer  program  is  thus  capable  of  treating  waveguides  that  are 
reasonable  model  representations  of  the  actual  complex  terrestrial  waveguide. 
The  program  does  not  allow  for  variation  of  the  guide  parameters  perpendicular 
to  the  propagation  path,  therefore  uniform  conditions  had  to  be  assumed 
over  distances  extending  a  Fresnel  Zone  to  either  side  of  the  path.  This 
restriction  has  long  been  recognized  as  a  major  limitation  of  the  model. 

2.2  Recent  Developments 

In  actual  propagation  cases  both  lateral  and  transverse  variations 
can  be  expected  to  occur  over  most  paths,  especially  at  high  latitudes. 
Recently,  off-path  effects  have  been  considered  and  important  steps  have 
been  taken  in  approaching  formulations  that  address  the  real  three-dimensional 


problems,  at  least  in  an  approximate  way  (Pappert,  1980;  Field  and  Joiner, 
1979;  Greifinger  and  Greifinger,  1978,  1979;  Booker,  1980)  as  briefly 
described  below. 

Pappert  (1980)  treating  localized  disturbances  in  a  surface  propagation 
model  coupled  with  a  moments  method,  finds  that  a  thick  ionization  patch 
at  E-region  heights  can  produce  6-8dB  fades  in  signal  strength  for  a  1.6Mm 
propagation  path.  Deepest  fades  occur  when  the  disturbed  area  is  well 
within  the  first  Fresnel  Zone,  and  a  varying  amount  of  fade  is  predicted 
depending  on  size  and  location  of  the  disturbance  area.  When  the  center  of 
the  disturbance  is  located  outside  the  second  Fresnel  Zone,  he  predicts 
that  most  of  the  amplitude  fades  are  damped  out. 

Field  and  Joiner  (1979)  use  an  integral-equation  approach  for  analyzing 
propagation  in  the  earth-ionosphere  waveguide  where  conditions  change  over 
distances  comparable  with  a  Fresnel  Zone.  They  derive  an  expression  for  the 
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relative  errors  introduced  by  neglecting  transverse  ionospheric  gradients 
over  the  path  and  find  that  full  wave  methods  must  be  applied  when  the 
effective  width  of  a  localized  disturbance  is  less  than  two-thirds  of  the 
width  of  the  first  Fresnel  Zone.  These  authors  also  conclude  that  the 
WKB  approximation  significantly  over-estimates  the  propagation  anomaly 
when  the  disturbance  is  centered  near  the  propagation  path  and  underestimates 
the  anomaly  when  the  disturbance  is  centered  far  off  path  (c.f.  also  Field, 
1978  ) • 

Further  physical  insights  into  ELF  propagation  and  ionospheric 
parameters  have  been  achieved  by  Greifinger  and  Greifinger  ( 1978, 

1979)  in  associating  the  eigenvalue  solution  mainly  with  ionospheric 
properties  at  two  well  defined  altitudes  where  maximum  ohmic  heating 
occurs.  Along  with  two  scale  heights  at  these  altitudes,  they  derive 
an  --’.pproximate  solution  for  the  ELF  mode  at  high  latitudes.  They  show 
that  the  lower  altitude  region  of  interest  occurs  where  the  conduction 
current  parallel  to  the  magnetic  field  becomes  equal  to  the  displacement 
current.  At  these  heights  the  associated  scale  height  is  that  of  the 
parallel  conductivity.  At  the  other  height  regions  the  scale  height  in 
question  is  that  of  the  refractive  index.  For  daytime  conditions  the 
authors  find  that  this  upper  height  level  is  situated  where  the  local 
wave  number  becomes  equal  to  the  reciprocal  of  the  local  scale  height  of  the 
refractive  index.  With  knowledge  of  these  four  parameters  an  approximate 
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wave  solution  may  be  found.  The  authors  conclude  that  the  principal  attenua¬ 
tion  mechanisms  are  Joule  heating  by  longitudinal  currents  in  the  vicinity 
of  the  lower  altitudes  and  energy  leakage  of  the  whistler  component  of  the 
ELF  wave  at  the  upper  altitudes. 

Booker  (1980)  presents  an  extensive  survey  of  ELF  propagation  and 
outlines  a  theory  for  propagation  in  the  earth-ionosphere  transmission  line 
based  on  the  work  of  Greifinger  and  Greifinger  (1979)  and  the  earlier 
theory  by  Booker  and  Lefeuvre  (1977).  His  method  greatly  advances  the  physical 
understanding  of  the  complex  reflection  and  absorption  processes.  Booker 
(1980)  draws  attention  to  several  results  important  for  understanding  ELF 
propagation,  among  others: 

-  At  ELF  the  transition  from  nearly  horizontal  phase  propagation  to 
nearly  vertical  propagation  takes  place  in  a  refracting  stratum  at  the 
bottom  of  the  ionosphere. 

Between  the  bottom  of  the  ionosphere  and  the  level  of  ELF  reflection, 
change  of  phase  can  as  a  first  approximation  be  calculated  on  the  basis 
of  vertical  propagation  in  free  space. 

-  The  reflecting  stratum  is  located  approximately  at  the  level  where  the 
ionosphere  fails  to  respond  as  a  doubly  refracting  medium  to  the  ELF 
waves . 

Thus,  above  the  reflecting  level  the  ionosphere  is  a  slowly-varying,  non¬ 
isotropic  medium,  whereas,  below  it,  the  ionosphere  behaves  essentially 
as  an  isotropic  medium. 

In  his  theory  Booker  includes  the  magnetic  field  and  derives  a  pro¬ 
pagation  model  faith  worldwide  application. 
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3.0  ELF  TRANSMITTER  AND  PROPAGATION  PATH  CHARACTERISTICS 

A  thorough  and  comprehensive  treatment  of  man-made  ELF  experiments 
is  given  in  a  special  issue  of  IEEE  Transaction  of  Communications  (April  1974). 
Limited  experimental  activity  has  continued  since  1974  using  the  US  Navy's 
experimental  transmitter  facility  in  Wisconsin. 

3.1  The  US  Navy  ELF  Transmitter  (Wisconsin  Test  Facility  -  WTF) 

The  Navy  ELF  Wisconsin  Test  Facility  (WTF)  is  located  in  the  Cheqeamegon 
National  Forest  in  north-central  Wisconsin,  USA,  and  consists  of  two  22.5  km 
north-south  antennas  (one  buried  and  one  elevated)  and  one  22.5  km  elevated 
east-west  antenna.  Each  antenna  is  grounded  at  both  ends.  The  transmission 
station  is  located  at  the  mid-point  intersection  of  the  two  antennas.  The 
antenna  array  pattern  can  be  steered  to  any  particular  receiving  direction. 

The  transmitter  has  been  operated  in  the  45  and  75  Hz  bands  with 
antenna  currents  of  300  A.  The  antenna  radiation  pattern  and  soil  con¬ 
ductivity  (ag)  under  the  antenna,  have  been  measured  by  Bannister  &  Williams 
(1974),  for  the  EW  antenna: 

—lx 

a  =  3.2  10  mho/m  at  75  Hz 
e 

-4 

a  =  2.8  10  mho/m  at  45  Hz 
e 

All  ELF  field  strength  values  presented  in  this  report  have  been 
normalized  to  the  above  antenna  characteristics,  and  the  revised  antenna 
radiation  pattern  (Bannister,  et  al.  1974)  is  applied. 

3.2  The  Ionospheric  Conditions  Near  the  Transmitter 

In  Figure  3.1  is  indicated  the  location  of  the  transmitter  as  well 
as  positions  of  selected  magnetic  L-shells.  The  Booker  excitation  circle 
(Booker,  1973)  is  plotted  for  a  frequency  of  75  Hz  (radius  equals  X^/211, 
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Figure  3.1  Location  of  the  Wisconsin  Test  Facility  (WTF)  and  positions  of 
the  auroral  zone  and  selected  magnetic  L-shells.  The  Booker 
excitation  region  is  also  shown  along  with  the  intensity  of 
precipitating  electrons  measured  on  a  satellite  pass.  See 
text  for  further  explanation. 
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where  Xq  is  the  free  space  ELF  wavelength).  Superposed  on  the  figure  is  the 
counting  rate  of  precipitating  energetic  electrons  (>  150  keV)  as  measured 
on  a  pass  by  the  1972-076B  satellite.  Such  latitude  variations  of  the  in¬ 
coming  fluxes  are  quite  typical,  and  will  cause  marked  local  perturbations 
in  the  ionosphere  below  100  km.  Thus,  in  this  case  homogeneous  ionospheric 
conditions  do  not  exist  over  the  whole  excitation  region  of  the  transmitter. 
Experience  from  several  years  of  satellite  measurements  indicate  that 
precipitating  electrons  in  the  L  range  of  3-5  are  probably  present  most  of 
the  time  exhibiting  dynamic  time  -  and  spatial  variations. 

Such  findings  complicate  the  ELF  field  prediction  schemes,  but  this 
problem  has  not  yet  been  attacked  fully  from  a  theoretical  point  of  view. 

In  model  calculations  homogeneous  conditions  over  the  transmitter  area  are 
assumed. 

3.3  ELF  Path  Characteristics 

The  NOSC  ELF  computer  code  allows  calculations  to  be  made  with  the 
propagation  path  divided  into  a  number  of  segments,  within  which  the  pro¬ 
pagation  medium  and  boundary  conditions  are  kept  constant.  Table  3.1  l^sts 
for  the  path  WTF  to  Tromso  the  segmentation  and  data  for  the  various  segments 
as  used  by  the  authors.  Other  paths  have  been  segmented  in  a  similar  manner 
(Imhof  et  al.,  1977). 

Effects  of  changes  in  the  ionospheric  conditions  due  to  different 
solar  zenith  angles  along  a  propagation  path  (local  time  dependence)  have 
also  been  considered  and  will  be  further  discussed  in  section  7. 
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3. 4  Ambient  Profiles  and  Collision  Frequencies 

For  the  ambient  day-  and  nighttime  ionsopheric  conditions  the  following 
electron  and  ion  density  profiles  have  been  used  by  the  present  authors. 

0  -  100  km  GE  Tempo,  RPT66TMP-82 

1  March  1966  by  W  .  Knapp 

100  -  200  km  Satellite  Environment  Handbook,  Stanford  University 
Press,  1965,  edited  by  F.  S.  Johnson 

These  ambient  profiles  are  reproduced  in  Figure  3.2  (from  Pappert  and 
Moler,  1974). 

The  electron-neutral  and  ion-neutral  collision  frequencies  used  both  for 
ambient  and  disturbed  conditions  are  likewise  those  listed  by  Pappert  and  Moler 
(1974),  cf  Table  3.2. 

Table  3.2  Collision  Frequencies 


Nighttime  Ionospheric  Collision  Frequency 
Altitude  Positive  Negative 

(km) _  Electrons _ Ions  _ Ions 


(sec- 

L) 

(sec'i) 

(sec'l) 

250 

I.05 

X 

102 

4.50 

x  10"1 

4.50 

x  10'1 

225 

3.50 

X 

10 1 

9.00 

x  10"1 

9.00 

x  10'1 

220 

3.00 

X 

10 1 

1.00 

x  10° 

1.00 

x  10° 

210 

3.30 

X 

10 1 

1.30 

x  10° 

1.30 

x  10° 

200 

4.50 

X 

10 1 

2.00 

x  10° 

2.00 

x  10° 

150 

1.60 

X 

10  3 

4.50 

x  101 

4.50 

x  101 

120 

1.00 

X 

104 

3.00 

x  102 

3.00 

x  102 

100 

3.90 

X 

io4 

8.00 

x  103 

8.00 

x  103 

0 

4.30 

X 

1011 

1.07 

x  10 10 

1.07 

x  10 10 
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While  these  electron  and  ion  density  profiles  have  been  used  in  previous 
work  on  this  project  and  throughout  the  present  work,  there  are  newer  daytime 
data  which  permit  a  more  suitable  daytime  profile  to  be  suggested  for  future 
work.  This  is  discussed  in  Section  10. 
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4.0  ELF  PROPAGATION  MEASUREMENTS  AND  PREDICTIONS 

A  number  of  reports  present  results  from  recordings  of  WTF  signals 
propagated  over  different  paths  in  the  USA  and  to  Europe  (Bannister  and 
Williams,  1974;  Bannister  1974;  Davis,  1976). 

Data  for  paths  such  as: 

-  WTF  to  Connecticut,  USA  (1.6  Mm  path  length) 

-  WTF  to  Greenland  (3.6  Mm  "  "  ) 

-  WTF  to  Tromso,  Norway  (6.0  Mm  "  "  ) 

-  WTF  to  Pisa,  Italy  (7.5  Mm  "  ”  ) 

have  been  collected  indicating  the  diversity  of  different  propagation 
media  analyzed. 

The  measurements  have  yielded  several  general  features  of  ELF 
propagation,  viz: 

-  Seasonal  changes  in  ionospheric  characteristics  influence  ELF 
propagation,  e.g.  through  variation  in  the  attenuation  rates  and 
relative  excitation  factors. 

Daytime  attenuation  rates  at  75  Hz  are  higher  than  nighttime  values, 
due  to  higher  electron  densities  below  80  km  altitude  resulting  from 

electron  detachment;  typical  ambient  values  are  1.1  -  1.3  dB/Mm  during 
the  day  and  0.8  -  1.0  dB/Mm  at  night. 

Nighttime  propagation  is  more  variable  than  daytime  propagation 

-  Polar  cap  and  auroral  zone  effects  are  found  in  ELF  propagation 

Such  findings,  reported  e.g.  by  Bannister  (1974)  are  also  reflected  in 
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theoretical  field  strength  predictions  using  propagation  models  with  ambient 
and  disturbed  ionospheres  (Pappert  and  Moler,  1974;  Pappert,  1980; 

Field  and  Joiner,  1979;  Imhof  et  al.,  1976,  1977,  1978,  1980). 


4.1  Short  Path  Propagation 

Initially  results  of  some  short  path  measurements  (<  2000  km  path 
length)  are  summarized.  Field  strengths  have  showed  that  anomalous  signal 
reductions  of  3  dB  and  more  may  occur  with  a  frequency  of  30%  of  all  nights. 

It  has  been  calculated  that  these  signal  strength  reductions  are  caused  by 
precipitation  of  electrons  from  the  radiation  belts  resulting  in  enhanced 
ionization  in  the  lower  atmosphere.  Attempts  to  correlate  the  anomalous 
signal  strength  behavior  with  geomagnetic  indicies  indicating  the  degree 
of  disturbed  magnetospheric  conditions  have  met  with  limited  success. 

As  an  example,  in  Figure  4.1  values  measured  by  Bannister  and  Williams 
(1974)  are  reproduced  over  the  years  1970-73  for  the  WTF  to  Connecticut  path. 
Seasonal  variations  are  seen,  but  the  average  value  of  all  observations 
agree  to  within  0.1  dB  of  the  predicted  ambient  night  value  as  calculated 
by  LMSC.  Furthermore,  in  Figure  4.2  is  depicted  an  anomalous  signal  reduction 
event  occuring  on  19  April  1973.  The  signal  reduction  lasting  several 
hours  amounted  to  :  5  dB  as  compared  with  night  values  before  and  after 
the  event.  Such  large  signal  attenuations  have  yet  to  be  recorded  at  a 
time  when  a  discrete  electron  precipitation  event  in  the  same  location  was 
measured  from  a  satellite. 
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Figure  4.2  16-20  April  1973  Connecticut  76-Hz  Nighttime  Field  Strengths 

versus  Local  Time  (Effective  Integration  Time  =  30  min/sample 
Bannister,  1973) . 
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4. 2  Long  Distance  Propagation 

Long  distance  ELF  propagation  (>1600  km)  also  exhibits  anomalous  signal 
fades.  Figures  4.3  and  4.4  show  measurements  over  long  and  short  paths  by 
Davis  (private  communication,  1977)  and  Bannister  (private  communication,  1978), 
respectively,  together  with  plots  of  the  magnetic  disturbance  index  Dst  and 
energetic  precipitating  electron  fluxes.  However,  no  simple,  systematic  rela¬ 
tionships  appear  to  exist  between  signal  fades  and  disturbed  conditions. 

4 . 3  Possible  Influence  of  Sporadic  E-Layers  on  ELF  Propagation 

Attempts,  however,  to  explain  the  signal  attenuation  due  to  excessive 
losses  in  sporadic  E-layers  have  resulted  in  predicted  signal  fades  of  the 
right  magnitude  (Pappert,  1978,  1980;  Barr,  1977;  Imhof  et  al.,  1976,  1977, 
1978,  1980)  both  for  short  and  long  paths. 

Measured  Eg  layers  are  very  thin  (3-4  km  or  thinner),  but  they  may 

2 

extend  over  large  geographical  areas,  up  to  1  Mm  and  more  (Whitehead,  1970). 
These  layers  do  not  reside  at  a  constant  height,  but  seem  to  move  towards 
lower  heights  with  time.  One  example  of  the  observed  shape  of  a  sporadic 
E  layer  measured  at  Wallops  Island  is  shown  in  Figure  4.3,  (Smith  and  Miller, 
1980).  As  noted  in  Appendix  A,  Eg  layers  are  observed  in  the  altitude  range 
100  -  130  km.  The  Eg  layer  has  been  placed  at  118  km  in  contrast  to  the 
observed  altitude  of  105  km  because  it  was  found  (Imhof,  et.  al . ,  1977)  that 
such  a  layer  produced  a  very  small  effect  at  105  km  and  a  large  effect  at  75  Hz 
at  118  km.  See  Appendix  A  for  further  general  characteristics  of  Eg  layers. 
Propagation  characteristics  at  75  Hz  for  such  a  profile  indicate 
significant  departures  from  ambient  conditions,  viz: 

Ambient  night  :  a  =  1.07  dB/Mm;  v/c  =  0.85 
Eg  :  a  =  1.83  dB/Mm;  v/c  =  0.75 
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eiving  station.  Also  plotted  are  the  fluxes  of  precipitating  electrons  >150  keV  measured 
the  S72-1  satellite.  The  D  geomagnetic  index  is  plotted  in  the  upper  section. 


Figure  4.4  Plot  of  the  precipitating  electron  flux  >160  keV  as  measured  from  the  satellite  on 
coordinated  passes  and  the  ELF  signal  strengths  as  received  at  the  various  stations. 
For  comparison  the  D  geomagnetic  index  is  also  plotted  to  indicate  the  general 
occurrence  of  magnetic  storms. 


MODIFIED  WALLOPS  ISLANDS  PROFILE 
Er  LAYER  MOVED  UP  TO  1 18  km 
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Figure  4.5  The  standard  nighttime  ambient  electron  density  profile  with  a  superimposed  experimental 
sporadic  E  layer  profile  taken  from  Smith  and  Miller  (1980),  but  placed  at  118  km. 
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For  a  thicker  ionization  ledge  positioned  at  123  km  altitude  Pappert 
(1980)  obtained  even  much  higher  attenuation  rates  for  a  short  path: 

a  =  9.83  dB/Mm;  v/c  =  0.68.  In  this  case  a  resonant  absorption  seems  to 
be  present  (cf  also  Imhof  et  al.,  1976,  1977;  Pappert,  1978).  However,  this 
resonance  is  strongly  dependent  upon  thickness  of  the  ionization  ledge  and 
its  position  in  altitude  (Imhof  et  al.,  1977).  Large  values  for  the  attenua¬ 
tion  rates  can  also  result  from  reflections  at  more  than  one  level  in  the 
atmosphere  (Barr,  1977). 

Amplitude  variations  also  on  longer  propagation  paths  may  tentatively 
be  explained  by  sporadic  E  layers,  in  as  much  as  propagation  characteristics 
for  such  profiles  indicate  attenuation  effects,  but  the  degree  of  attenuation 
depend  on  the  position  of  the  Eg  layer  with  respect  to  the  propagation  path. 
As  an  example,  in  Table  4.1  is  shown  some  computed  field  strength  values  for 
Eg  conditions  over  three  selected  areas  for  the  path  WTF  to  Tromso,  Norway. 
The  values  marked  Eg  x  10  are  calculated  for  an  increase  in  the  electron 
density  of  the  Eg  layer  (only)  by  a  factor  of  10.  Apparently  the  geographic 
position  of  the  Eg  layer  with  respect  to  the  propagation  path  and  its 
intensity  must  be  known  in  order  to  predict  ELF  field  strength  values  for 
comparison  with  measured  values. 
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Table  4.1  Calculated  field  Strengths  at  75  Hz  for  propagation  WTF  to  Tromso. 


Field  Strength  (dB  WRT  lA/m)  * 


Ambient  Night  -154.7 
Eg  over  receiver,  ambient  elsewhere  -155.5 
Eg  over  Norwegian  Sen,  ambient  elsewhere  -156.5 
Eg  x  10  over  Norwegian  Sea,  ambient  elsewhere  -165.9 


*  The  field  strength  values  in  this  table  and  in  all  other  later  LMSC 

references  in  this  report  have  been  normalized  to  include  the  signal  of  the 

WTF  E-W  antenna  in  addition  to  that  of  the  N-S  antenna  which  is  computed 

by  the  waveguide  mode  code:  The  E-W  antenna  conditions  are:  Current  300A; 

-4 

length  22.5  km,  conductivity  of  ground  under  antenna  3.2  x  10  mho/m; 
phasing  between  E-W  and  N-S  antennas  =  0°. 

Table  4.2  shows  in  more  detail  examples  of  computed  characteristics 
at  75  Hz  for  propagation  from  WTF  to  Tromso  for  ambient  and  sporadic  E 
conditions. 

Effects  are  predicted  for  an  E  over  Norwegian  Sea  area  ('v  1400  km 

in  lateral  size)  amounting  to  2  dB  attenuation.  For  the  very  enhanced 

6  3 

Eg  (Peak  electron  density  of  1.5  .  10  el/cm  )  a  major  signal  fade  in 
excess  of  10  dB  is  predicted. 

We  may  thus  conclude  that  idealized  Eg  disturbances  have  produced 
severe  ELF  attenuation  in  model  calculations,  but  actual  measurements  of 
Eg  conditions  have  not  been  made  during  the  operation  of  the  Navy's  ELF 
transmitter.  Attempts  to  explain  the  observed  ELF  signal  fades  in  terms 
of  absorption  due  to  Eg  conditions,  can  therefore  not  be  conclusive,  but  the 
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Table  4.2  Propagation  Characteristics  at  75  Hz  calculated  for  ambient  and 
sporadic  E  conditions  for  the  WTF  area  (o  *  3. 1C  ~4  mho/m) , 
Greenland  Ice  caps  (a  -  1.10”^  mho/m)  and  Norwegian  Sea 
(a  -  4  oho/m). 


I 

I 

I 

I 


Conditions 

Eigenangle  (Degrees) 

S  -  sin  6 

a 

v/c 

Ambient 

(day) 

(81.918  -  1  38.061) 

(1.22  -  1  0.10) 

1.37 

0.82 

Ambient 

(night) 

(82.795  -  1  34.328) 

(1.18  -  1  0.08) 

1.09 

0.85 

Es 

(O-3.10-4) 

(82.158  -  1  46.395) 

(1.33  -  1  0.12) 

1.68 

0.75 

Es 

(a-1.10-5) 

(79.549  -1  49.804) 

1.38  -1  0.18) 

2.43 

0.73 

Es 

(o-4) 

(80.191  -1  46.935) 

1.33  -1  016) 

2.13 

0.75 

Es  x  10 
(a-3.10"4) 

(78.041  -1  58.076) 

(1.52  -1  0.25) 

3.39 

0.66 

Es  x  10 
(o-1.10~5) 

(69.888  -1  62.374) 

(1.55  -  1  0.45) 

6.18 

0.64 

Es  x  10 
(o-4) 

(59.926  -  1  58.767) 

(1.36  -  1  0.61) 

8.54 

0.73 
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theoretical  efforts  in  this  area  point  out  the  potential  influences  of  Eg 
on  an  ELF  communication  system. 

4.4  Propagation  During  REP  Conditions 

Some  electron  precipitation  events  are  classified  as  Relativistic 
Electron  Precipitation  (REP)  events  since  it  is  believed  that  electrons  with 
energies  between  0.1  and  1  Mev  are  principally  responsible.  These  events 
cause  significant  ionization  in  the  lower  D  region  below  about  70  km.  The 
events  are  observed  to  be  strongest  in  the  region  below  the  edge  of  the 
outer  radiation  belt  during  the  daytime  hours  and  are  believed  to  result 
from  the  scattering  of  trapped  relativistic  electrons  by  some  particle 
interaction.  Substorm  activity  is  believed  to  be  a  necessary  condition  for 
REP  events  but  not  all  substorms  lead  to  internal  relativistic  electron 
precipitation. 

Relativistic  electron  precipitation  events  frequently  cause  excess 
ionization  in  the  ionospheric  D-region  to  persist  for  one  to  several  hours 
over  relatively  large  geographical  areas.  Larsen  (1974,  1979)  infers 
that  the  disturbed  area  may  extend  over  2000  km  in  longitude.  Satellite 
measurements  indicate  that  the  precipitation  area  may  cover  several  degrees 
in  latitude  and  one  may  therefore  assume  that  the  REPs  may  have  such  large 
geographic  extentions  that  they  can  influence  ELF  propagation. 

Examples  of  precipitating  electron  energy  spectra  are  shown  in  Figure 
4.6.  Note  the  "peaks"  at  energies  between  800-1500  keV.  Ion  pair  production 
rates  from  such  events,  are  high  enough  to  give  substantial  increases 
in  the  electron  density  down  to  50-60  km  in  altitude,  cf  Figure  4.7.  Assuming 
that  this  precipitation  occurred  over  segments  of  the  propagation  paths  under 
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study,  as  indicated  in  Figure  4.8,  it  is  possible  to  estimate  the  effects 
on  propagation,  as  illustrated  in  Tables  4.3  and  4.4.  It  is  seen  that  signal 
enhancements  are  predicted  for  all  these  nighttime  REPs,  by  1  dB  or  more  when 
the  effects  occur  over  the  transmitter  (receiver)  or  tenths  of  dB  when  the 
effects  are  localized  to  areas  over  the  paths.  REP  conditions  over  the  com¬ 
plete  path  would  markedly  enhance  the  signal  (by  several  dB) . 

Propagation  characteristics  calculated  at  75  Hz  for  ambient  conditions 
and  for  two  of  the  REP  events  of  Figure  4.7,  in  three  segments  of  the  WTF- 
Tromso  path,  are  shown  in  Table  4.5. 

In  some  cases  REP  ionization  profiles  have  also  resulted  in  predicted 
signal  decreases  depending  upon  the  relative  position  between  the  disturbed 
area  and  the  propagation  path  (Imhof  et  al .  1978). 

In  addition  to  direct  ionization  due  to  the  electrons  during  an  REP 
event,  there  is  also  some  secondary  ionization  due  to  bremsstrahlung  x-rays 
produced  by  the  electrons.  This  contribution  is  usually  neglected,  but  in 
some  cases  it  should  be  considered.  In  Figure  4.9,  it  is  shown  that  the 
x-rays  may  dominate  as  an  ionization  source  between  40  and  60  km,  filling  in 
the  "valley  between  the  profiles  of  the  ionization  by  electrons  and  that 
by  galactic  cosmic  rays.  Inclusion  of  the  x--ray  contribution  resulted 
in  a  predicted  field  strength  attenuated  by  0.7  dB  over  the  signal  strength 
computed  without  the  x-ray  ionization  (calculated  for  the  1.6  Mm  path:  VJTF  to 
Connecticut) . 

In  summarizing,  ELF  measurements  coordinated  with  satellite  particle 
measurements  and  ELF  computer  code  predictions,  have  given  the  following 
conclusions : 
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Table  4.3  Calculated  Field  Stre  ?ths  for  Electron  Precipitation  Events. 
Transmissions  from  the  Wisconsin  Test  Facility  to  the 
Receiving  Station  at  Tromso,  Norway* 


Event  Date 

Observed  Invariant 
Latitude  Range  of 
Precipitation 

Path  Segment 
for  Precipitation 

Calculated  Field  Strength 
at  75  Hz  WTF  -  Tromso 
(dB  wrt  to  lA/m) 

Ambient  Night 

- 

- 

-  154.7 

26  Mar.  1976 

57.7°  -  64.8° 

0  -  900  km 

-  153.6 

19  Apr.  1977 

61.5°  -  65.7° 

0  -  900 

-  153.5 

200  -  900 

-  154.6 

whole  path 

-  151.9 

7  May  1977 

65.2°  -  68.9° 

0  -  1800 

-  153.4 

1000  -  1800 

-  154.4 

whole  path 

-  151.9 

*  The  whole  oath  length  is  6000  km 


Table  4.4  Calculated  Field  Strengths  WTF  -  Greenland*  for  19  April  and  7  May  1977 
Precipitating  Electron  Profiles 


Profile 

Date 

Path  Segments 
for  Precipitation 

Calculated  Field 
Strength  at  75  Hz 
(dB  wrt  1  A/m) 

Calculated  Change 
from  Ambient  Night 
(dB  wrt  1  A/m) 

19  April 

0  -  1200  km 

-  152.1 

+  1.2 

300  -  1200  km 

-  152.2 

+  0.1 

7  May 

0  -  1200  km 

-  151.2 

+  1.1 

300  -  1200  km 

-  152.1 

+  0.2 

Ambient 

Night 

-  152.3 

*  Total  path  length  3500  km 
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Table  4.5  Propagation  characteristics  at  75  Hz  calculated  for  ambient  and  REP 
conditions  for  the  Wit  area  (o  =  3.10  4  ntho/rn),  Greenland  icecap 
(o  -  1.10  mho/m)  and  Norwegian  Sea  (a  =  4  mho/m). 


Conditions 

Eigenangle 

(9),  Degrees 

S  =  sin  9 

a 

(db/Mm) 

v/c 

Ambient  (day) 

(81.918  - 

i 

38.061) 

(1.22  -  i 

0.10) 

1.37 

0.82 

Ambient  (night) 

(82.795  - 

i 

34.328) 

(1.18  -  i 

0.08) 

1.09 

0.85 

REP  26  Mar  76 
(a-3.10-4) 

(83.772  - 

i 

33.313) 

(1.17  -  i 

0.07) 

0.91 

0.86 

REP 

(a=1.10“5) 

(80.646  - 

i 

38.667) 

(1.22  -  i 

0.12) 

1.61 

0.82 

REP 

(o=4) 

(84.729  - 

i 

31.896) 

(1.15  -  i 

0.05) 

0.74 

0.87 

REP  19  Apr  77 
(a=3. 10-4) 

(83.116  - 

i 

33.036) 

(1.16  -  i 

0.07) 

1.00 

0.86 

REP  . 
(a-1.10“5) 

(80.027  - 

i 

33.036) 

(1.22  -  i 

0.13) 

1.73 

0.82 

REP 

(a=4) 

(84.054  - 

1 

31.401) 

(1.15  -  i 

0.06) 

0.81 

0.87 
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From  both  the  satellite  and  ELF  station  measurements  it  has 
been  found  that  direct  particle  precipitation  into  the  atmosphere 
can  cause  ELF  transmission  anomalies.  In  these  anomalies  the  signal 
strengths  may  be  either  attenuated  or  enhanced  depending  upon  the 
geometry  and  details  of  the  ion-pair  production  profiles.  Nighttime 
REP  events  may  on  the  average  improve  ELF  propagation  conditions. 
Variations  in  the  nighttime  ELF  signal  strengths  on  a  fine  time  scale 
are  observed  which  may  be  due  entirely  to  electron  precipitation,  but 
cannot  be  accounted  for  quantitatively  due  to  present  limitations  in 
the  measurements  and  computational  techniques. 

The  geometry  for  the  effect  of  electron  precipitation  on  nighttime 
ELF  transmission  is  very  complex  and  as  a  result  the  need  for  new 
techniques  for  mapping  electron  precipitation  profiles  simultaneously 
over  a  broad  spatial  region  has  become  apparent. 
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5.0  CALCULATIONS  FOR  SIMULATED  PCA  CONDITIONS 

There  have  been  no  actual  measurements  of  signals  from  WTF  during 
a  PCA  event.  Calculations  for  simulated  conditions,  however,  indicate 
that  significant  signal  strength  reductions  could  occur.  The  total  amount 
of  signal  absorption  would  depend  on  the  extent  of  the  disturbed  area 
and  location  of  the  ELF  transmitter  with  respect  to  the  PCA  area.  Imhof 
et  al.  (1980)  summarized  their  findings  as  follows: 

-  Typical  PCA  nighttime  conditions  as  observed  during  the  August 
1972  event,  will  decrease  the  ELF  signal  at  75  Hz  by  several  dB 
over  a  polar  propagation  path  (WTF  to  Tromso)  for  uniform  conditions 
along  the  whole  path. 

-  The  received  ELF  field  strength  for  such  a  path  is  predicted  to  be 
reduced  by  an  additional  1-3  dB  when  the  WTF  transmitter  area  is 
outside  the  PCA  disturbed  area. 

The  ELF  signal  strength  was  found  to  be  sensitive  to  changes  in  the 
ion  densities  at  heights  below  45-50  km. 

5 . 1  PCA  Simulations  Using  Data  from  the  August  1972  PCA 

In  the  following  we  present  some  additional  work  (including  daytime 
PCA  simulations)  on  signal  strength  calculations  during  simulated  PCA 
conditions.  For  reference  the  PCA  profiles  used*  are  reproduced  in  Figures 

*  These  test  cases  are  based  on  electron  density  measurements  made  by  the  SRI 
Chatanika  incoherent  scatter  radar.  As  originally  interpreted  the 
measurements  yielded  electron  densities  which  are  now  known  to  be  too 
large  by  a  factor  which  is  altitude  dependent  and  ranges  in  magnitude 
from  one  to  two.  Theory  developed  by  Matthews  (1978)  and  Fukuyama 
and  Kofman  (1979)  shows  that  when  negative  ions  are  taken  into  account 
smaller  electron  densities  result  from  the  measurements  as  discussed 
in  Imhof  et  al.  (1980)  The  positive  ion  densities  inferred  from  the  electron 
densities  are  thus  too  small.  These  test  cases  however  have  been  used 
frequently  in  our  work  and  are  retained  as  good  examples  of  a  very 
intense  PCA. 
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5.1  and  5.2  (cf  Imhof  et  al.,1980  for  further  comments).  Care  has  been 
taken  to  include  ion  densities  at  all  heights  commensurate  with  the  electron 
densities . 

The  values  for  the  positive  ions  above  40  km  were  obtained  using  the 
chemistry  code  of  Gunton  et  al.  (1977).  At  lower  altitudes  the  values  for 
the  positive  ion  density  are  obtained  using  the  formula 

N+  =  ✓q/ai(h,T) 

where  height  and  temperature  dependent  values  of  the  ion-ion  recombination 

rate  a.  are  used,  of  Table  5.1 
x 


Table  i.l  Ion-ion  Recombination  Rates  (a)  Calculated 
for  65°N,  August  ' 


Altitude 

(km) 

P 

Pressure 

(Dynes/cm2) 

T 

Temperature 

<K) 

3  - 1 

(cm  sec  ) 

60 

2.869  102 

261 

5.16  10“8 

55 

5.366  102 

274 

5.26  10"8 

50 

9.834  102 

278 

5.46  10~® 

45 

1.801  103 

274 

5.87  10~8  1 

40 

3.383  103 

261 

6.84  10"8 

35 

6.571  103 

247 

9.11  10"8 

30 

1.305  104 

236 

1.41  10~7 

25 

2.800  104 

228 

2.64  10~7 

20 

5.900  104 

225 

5.16  10“7 

15 

1.260  105 

225 

1.05  10'6 

10 

2.680  105 

225 

2.17  10'6 

5 

5.410  105 

260 

3.03  10'6 

0 

1.017  106 

289 

4.35  10‘6 

Note:  These  values  for  ^  have  been  calculated  using  the  formula 
=  5  x  10'8  +  6  x  10'6  P  T'5'/2, 

where  P  is  pressure  in  dynes/cm2  and  T  is  temperature  in  degrees  Kelvin. 

This^ formula  is  from  Niles  (1976)  with  the  first  term  changed  from  3  x  10‘8  to 
5  x  10  in  light  of  the  latest  laboratory  measurements  by  Smith  and  Church  (1977). 
the  model  atmosphere  is  for  August  at  65°N  latitude. 
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Table  5.2  Comparison  between  ELF  absolute  field  strength  calculations  at  75 
Hz  for  the  path  WTF  to  Tromso  . 


LMSC 

Predicted 

Condition 

dB/lA/m 

Change  (dB) 

Day,  ambient 

-154.7 

- 

Night,  ambient 

-154.7 

" 

Day  PCA  Hi 

-160.4 

-5.7 

Day  PCA  HI 

-157.5 

-2.8 

Day  PCA  #3 

-157.0 

-2.3 

Night  PCA  HI 

-161.1 

-6.4 

Night  PCA  HI 

-162.6 

-7.9 

5.2  Sensitivity  of  ELF  Field  Strengths  Using  Other  Simulated  PCA  Conditions 
The  earlier  efforts  (Imhof  et  al.,  1976,  1977,  1978,  1980)  have 
emphasized  the  importance  of  the  conditions  prevailing  during  the  August 
1972  event  in  absorbing  energy  from  the  ELF  waves.  The  present  work  extends 
these  efforts  through  use  of  data  for  three  PC...  i  vents,  having  energy  spectra 
of  different  spectral  shapes  and  intensities,  viz: 

PCA  of  12  Nov  1960  1930  UT 

PCA  of  2  Nov  1969  1448  CST 

PCA  of  4  Aug  1972  Daytime  Test  Case  #1 

The  ion  and  electron  density  profiles  for  the  PCA  of  12  November  1960 
were  calculated  using  the  ion  pair  production  altitude  profiles  shown  in 
Figure  5.3  together  with  electron  and  ion  loss  coefficients  calculated 
during  our  modeling  of  the  PCA  of  4  August  1972  (Imhof  et  al.  1980). 
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The  ion  and  electron  density  profiles  for  the  PCA  of  2  November  1969 
were  calculated  with  our  ion  chemistry  model.  The  ion  pair  production  alti¬ 
tude  profile  time  dependence  was  derived  from  Swider  (1975) . 

The  deduced  electron  and  positive  ion  densities  for  these  cases  are  given 

in  Table  5.3  and  plotted  in  Figure  5.4.  It  may  be  noted  in  Figure  5.3 

that  the  November  1960  event  had  a  harder  proton  spectrum  yielding  higher 

N+  densities  below  30  km  than  the  August  1972  PCA.  Above  80  km,  however, 

the  ionization  due  to  the  protons  during  the  1960  and  1969  events  was  below 

the  normal  ionization  levels  of  the  solar  radiation,  and  the  daytime  ambient 

N  and  N+  curves  have  been  used  from  80  to  160  km  for  these  two  cases, 
e 

We  first  present  the  calculated  field  strengths  at  75  Hz  for  propagation 
from  WTF  to  Tromso  which  are  as  follows: 

PCA  12  Nov  1960  -163.3  dB;  8.6  dB  below  ambient 

PCA  2  Nov  1969  -155.7  dB;  1.0  dB  below  ambient 

PCA  4  Aug  1972  -160.4  dB;  5.7  dB  below  ambient 

(Day  test  case  It 2) 

For  these  simulated  cases,  ambient  conditions  were  assumed  over  the 
WTF  transmitter  area.  In  the  calculation  for  the  PCA  of  2  November  1969, 
the  standard  segmentation  of  the  path  into  eight  parts  was  used.  The  first 
five  segments  were  assumed  to  be  in  daylight  and  the  remainder  in  night. 

For  the  daylight  segmet.es  beyond  the  WTF  area,  a  daylight  profile  calculated 
for  segment  3  under  PCA  conditions  was  used.  The  segments  assumed  to  be 
in  darkness  were  given  a  profile  based  on  nighttime  PCA  calculations  for 
segment  3.  In  the  other  two  cases  all  of  the  segments  beyond  the  WTF  area 
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Table  5.3  Electron  and  Ion  Densities  for  Selected  PCA  Events 


Altitude 

(k.n) 

PCA  12 

Nov  60 

PCA  2 

Nov  69 

PCA  4 

Aug  72 

Ne(cm-~* ) 

N+(cm  3 ) 

Ne(cm“3) 

N+(cm“3) 

Ne(cm-3) 

N+(cm-3) 

100 

6. 00 (+4 ) 

6. 00 (+4) 

6.00(4-4) 

6.00(4-4) 

5 . 00 (+5 ) 

5. 00 (+5) 

95 

3.60(+4) 

3.60(+4) 

3.60(4-4) 

3.60+4) 

1.00(+5) 

1.00(+5) 

90 

1.50 (+4) 

1.50(+4) 

2.00(4-4) 

2.00(+4) 

6. 50 (+4) 

6. 50 (+4) 

85 

7. 00 (+3) 

7.00(4-3) 

1.57(4-4) 

1.57  (4-4) 

4. 80 (+4) 

4.80(+4) 

80 

5. 36(+3 ) 

5.60(4-3) 

1.74(4-4) 

1.75 (+4 ) 

3.80(4-4) 

3. 80 (+4) 

75 

6.40(+3) 

8.19(4-3) 

1.06(4-4) 

1.09(4-4) 

5. 00 (+4) 

5. 00 (+4) 

70 

6.46(+3) 

1.40(4-4) 

1.01(4-4) 

1.23 (+4) 

5. 50 (+4) 

7.00(4-4) 

65 

5.72(+3) 

3.04(4-4) 

9.22(4-3) 

2.07 (+4) 

5. 40 (+4) 

1.50(+5) 

60 

3.84(+3) 

6.01(4-4) 

7.32(4-3) 

3.17 (+4 ) 

4. 50 (+4) 

2.50(+5) 

55 

2. 18(+3) 

1.06(4-5) 

3.81(4-3) 

6.38(4-4) 

2.30(+4) 

3. 70(+5) 

50 

8. 50 (+2 ) 

1.59(4-5) 

1.68(4-3) 

9.18(+4) 

1.40 (+4 ) 

5.00(+5) 

45 

5. 46 (+2) 

2.08(4-5) 

2.39(4-2) 

8.84(+4) 

1.10(+3) 

6. 00 (+5) 

40 

l.ll(+2) 

2.43(4-5) 

2.18(4-1) 

5. 69 (+4) 

2. 15 (+2) 

6. 00 (+5) 

35 

2.10(+1) 

2.47(4-5) 

9. 20 (-4) 

3.25(+3) 

1.00(+1) 

5.40(+5) 

30 

1.10(+1) 

2.22(4-5) 

3.00(-4) 

3. 55(+3) 

1.67 (-4) 

3. 75(4-5 ) 

25 

5 . 00 ( — 3 ) 

1.76(4-5) 

1.00 (-4) 

3.69(+3) 

7 . 70(-5) 

1.15(+5) 

20 

1.00 (-3) 

1.19(4-5) 

3.43(-5) 

3. 65(+3) 

3. 43 (-5) 

2. 05 (+4) 

15 

3. 00 (-4) 

6.14(4-4) 

1 . 50(-5) 

3.63(+3) 

1.50(-5) 

3. 63 (+3) 

10 

6. 58(-6 ) 

3.72(4-3) 

6.58(-6) 

3.72(+3) 

6. 58(-6) 

3.72(4-3) 

5 

3.43(-6) 

4.48(4-3) 

3.43(-6) 

4. 4 8 (+3) 

3.43(-6) 

4.48(+3) 

NOTE:  In  the  calculations  the  profiles  for  12  Nov  60  and  2  Nov  69  were 

extended  to  160  km  with  ambient  daytime  values  above  80km  and  100  km,  respec¬ 
tively.  The  profile  for  PCA  2  Nov  69  is  for  daytime  conditions. 
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used  the  profile  specified  in  Figure  5.3.  In  Table  5. A  are  listed  the  S-values 
(sin  0,  the  eigenangle  solution)  for  three  selected  segments  of  the  path, 
representing  a  range  in  ground  conductivity  and  different  propagations  with 
respect  to  the  earth's  geomagnetic  field.  For  comparison  a  recent  S-value 
computed  by  Field  et  al.  (1980)  at  75  Hz  for  a  "strong  PCA"  case  is  included. 

Detailed  information  on  the  computed  attenuation  rates  and  phase 
velocities  are  given  in  Table  5.5. 

Figure  5.5  shows  the  variation  of  ELF  field  strengths  as  a  function 
of  distance  for  these  three  cases. 


Table  5.4  Computed  values  for  S  *  sin  0  at  75  Hz  for  various  segments  along 

propagation  path  WTF  to  Tromso.  Also  given  is  the  value  deduced  bv 
Field  et  al.  (1980)  for  a  "strong  PCA"  case. 


Segment  PCA  12  Nov  60  PCA  2  Nov  69  PCA  4  Aug  72 

(Day  test  case  # 1) 

//6  Greenland  Icecap  1.4526  -  0.3084  i  1.3563  -  0.1916  i  1.4308  -  0.2725  i 

o  =  10  ^  mho/m 

# 7  Norwegian  Sea  1.3656  -  0.2139  i  1.2681  -  0.1104  i  1.3337  -  0.1725  i 

a  =  4  mho/m 

#8  Tromso  Area  1.3748  =  0.2252  i  1.2767  -  0.1187  i  1.3433  -  0.1830  i 

a  =  10  mho/ra 


Field  et  al  (1980)  Sj  -  1.50  -  0.25  i 

Value  for  "strong  PCA" 
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Table  5.5 

Propagation  characteristics  at  75 
different  proton  energy  spectral 

Hz  for  three 
shape. 

PCA  events 

of 

Segment  // 

PCA 

Nov  60 

PCA 

Nov 

69 

PCA  Aug  72 

a 

v/ c 

a 

v/c 

a 

v/c 

(db/Mm) 

(db/Nta) 

(db/ftn) 

1 

2.614 

0.7400 

2 

2.940 

0.7321 

1.641 

0.7828 

2.496 

0.7445 

3 

2.811 

0.7369 

1.527 

0.7882 

2.354 

0.7499 

4 

2.935 

0.7318 

1.637 

0.7829 

2.496 

0.7444 

5 

2.801 

0.7363 

1.520 

0.7884 

2.353 

0.7499 

6 

4.032 

0.6912 

2.612 

0.7373 

3.715 

0.6989 

7 

2.782 

0.7356 

1.505 

0.7886 

2.352 

0.7498 

8 

2.912 

0.7306 

1.618 

0.7833 

2.495 

0.7444 

5.3  Sensitivity  of  ELF 

Waves  to  Absorption 

Ions  at  Low 

Altitudes 

In  this  section  we  will  briefly  discuss  the  importance  of  ions  below 
45-50  km  altitude  in  determining  the  absorption  of  ELF  waves.  Calculations 
were  made  at  75  Hz  for  the  WTF  to  Tromso  path.  Ambient  day  conditions  were 
assumed  over  the  WTF  area. 

To  achieve  this  we  have  modified  the  PCA  4  Aug  1972  and  12  Nov  1960 
ion  profiles  at  low  altitudes,  and  calculated  the  effects.  Figures  5.6 
and  5. 7, respectively,  show  the  modifications  and  the  corresponding  field 
values . 

It  can  be  concluded  that  the  ions  formed  by  the  high 
energy  tail  of  the  precipitating  protons  below  45  -  50  km  are  significant. 
In  the  case  of  PCA  Nov  1960  calculation  shows  that  the  ions  below  40-45  kir. 
would  cause  about  8  dB  attenuation  over  the  6  megameter  path  to  Tromso, 
a  very  large  effect.  The  region  20  -  45  km  seems  to  be  the  most 
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Figure  5.7  Electron  and  ion  density  profiles  for  the  PCA  of  12  November  1960  along  with  modifications 

in  the  positive  ion  portions  to  test  sensitivity  of  the  signal  strength  to  such  modifications. 
The  numbers  on  the  profiles  are  attenuations  in  dB  at  75  Hz  over  the  WTF  to  Tromso  path. 
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important  height  range  for  this  attenuation,  but  enhanced  densities  even 
below  20  km  cause  some  absorption 

Figure  5.6  showing  the  modifications  to  the  daytime  4  Aug  1972  PCA 
(Test  Case  #1),  also  indicates  effects  due  to  a  postulated  "extreme"  case 
for  ion  densities  below  20  km.  Here  the  ion  densities  are  almost  two  orders 
of  magnitude  above  ambient  at  the  20  km  level.  The  amount  of  ohmic  heating 
that  occurs  at  these  altitudes  is  strongly  reduced  and  only  a  minor  ('v  ldB) 
amount  of  excess  absorption  takes  place.  Detailed  propagation  characteristics 
are  listed  in  Tables  5.6  and  5.7. 

From  the  present  and  earlier  studies  the  following  conclusions  can 
be  drawn  regarding  importance  of  ions: 

-  Relatively  little  ELF  attenuation  seems  to  occur  below  about 

20  km.  For  extreme  cases  with  high  N+  densities  some  ELF  attenuation 
is  predicted  to  occur  even  at  these  altitudes. 


Table  5.6  Propagation  Characteristics  at  75  Hz  for  PCA  day  Test  Case  It  1 
(A  Aug  72)  -  Original  and  Modifications. 


Segment 

it 

Original 

Profile 

Original  to 

30  km;  below 

30  km:  PCA 

Nov  60 

Original  to 

30  km;  Very  High 
Ion  Densities 

Original  to 

45  km;  Below 

45  km:  Sharp 
Transit . 

a(dB/Mkn) 

v/c 

a 

v/c 

a 

v/c 

a 

v/c 

1 

2.61A 

0.7A00 

2 

2.A96 

0.7AA5 

2.6AA 

0.7A52 

2.869 

0.7460 

1.141 

0.8036 

3 

2.35A 

0.7A99 

2.501 

0.7505 

2.724 

0.7513 

1.018 

0.8099 

A 

2.A86 

0.7AAA 

2.6AA 

0.7A52 

2.869 

0.7460 

1.141 

0.8036 

5 

2.353 

0.7A99 

2.500 

0.7505 

2.723 

0.7513 

1.018 

0.8098 

6 

3.715 

0.6989 

3.873 

0.6999 

4.113 

0.7011 

2.185 

0.7516 

7 

2.352 

0.7A98 

2.A99 

0.7505 

2.722 

0.7513 

1.017 

0.8098 

8 

2.A95 

0.7AAA 

2.643 

0.7451 

2.868 

0.7460 

1.140 

0.8036 
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Table  5.7  Propagation  characteristics  at  75  Hz  for  PCA  12  Nov  60  profile 
(original  and  modifications). 


Segment 

# 

Original 

Profile 

Original 
below  30 
day  test 

to  30  km 
km :  PCA 
case  ill 

Original  to  45  km 
below  45  km:  Sharp 
transition 

a 

(db/Mm) 

v/c 

a 

(db/Mm) 

v/ c 

a 

(db/Mm) 

v/c 

1 

3.039 

0.7278 

2.908 

0.7269 

2 

2.940 

0.7321 

2.810 

0.7312 

1.795 

0.7584 

3 

2.811 

0.7369 

2.682 

0.7360 

1.678 

0.7656 

4 

2.935 

0.7318 

2.804 

0.7309 

1.789 

0.7581 

5 

2.861 

0.7363 

2.672 

0.7354 

1.668 

0.7630 

6 

4.032 

0.6912 

3.895 

0.6901 

2.780 

0.7142 

7 

2.782 

0.7356 

2.653 

0.7347 

1.648 

0.7624 

8 

2.912 

0.7306 

2.782 

0.7297 

1.766 

0.7569 

-  In 

the  height 

range  20- 

■45km  ions  dominate  as 

the  cause  o 

f  ELF  wave 

attenuation  through  ohmic  dissipation  during  disturbed  conditions. 

Such  ohmic  losses  may  account  for  50-60%  of  the  total  wave  absorption, 
and  they  maximize  around  30-40  km  for  the  profiles  analyzed. 

At  heights  above  45-50  km  the  ion  contribution  to  the  ionospheric 
conductivity  decreases  fast  relative  to  the  electron  component 
which  rapidly  becomes  the  dominant  part  at  these  and  higher  altitudes. 
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5.4  Sensitivity  of  ELF  Waves  to  Changes  in  Electron  Densities  at  Low 

Heights  During  PCA  Conditions 

During  PCA  daytime  conditions  the  electron  densities  around  40  to  50  km 
may  be  4-5  orders  of  magnitude  above  the  ambient  values.  Clearly  some 
uncertainty  will  be  connected  with  the  deduced  PCA  densities  due  to  changes 
in  reaction  rate  coefficients  .  A  simple  test  on  the  sensitivity  of 
the  ELF  field  strengths  to  Ng  changes  at  low  heights  can  be  made  by  moving 
the  lower  part  of  the  Ng  profile  up  and  down,  say  5  km,  and 
strength  calculations  for  these  cases. 

As  an  example  we  have  taken  the  PCA  Day  Test  case  #1  and  modified  it 
below  60  km  in  the  way  suggested  above.  The  N+  curve  was  kept  unchanged 
as  were  the  electron  densities  above  65-70  km  (cf  Figure  5.8). 

Results  for  propagation  at  75  Hz  from  WTF  to  Tromso  were  as  follows: 

PCA  Daytest  case  //I 

Unchanged  Up  5  km  Down  5  km 

-  160.4  dB  -  161.3  dB  -  159.4  dB 

These  calculations  indicate  that  for  a  solar  proton  event  like  that 
of  Aug  72  an  uncertainty  in  the  lower  part  of  the  Ng  profile  of  +  5  km, 
would  yield  +  1  dB  uncertainty  in  the  field  strength  value.  In  a  real 
situation  the  Ng  values  would  probably  be  known  to  a  greater  accuracy  than 
assumed  in  this  example . 

It  may  be  noted  that  lowering  the  profile  at  low  altitudes  by  5  km 
(making  the  ionosphere  5  km  thicker)  increases  the  expected  signal  level 
essentially  through  a  lower  value  for  the  attenuation  constant.  Conversely, 
raising  the  Ng  profile  by  5  km,  resulted  in  increased  attenuation  rates  by 
0. 3-0.4  dB/Jfa. 
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5.5  Criteria  for  ELF  Field  Program 

Earlier  results  (Imhof  et  al. ,  1980)  indicated  that  the  predicted 

field  strength  at  75  Hz  for  WTF  to  Tromso  decreased  in  approximately  a 

logarithmic  manner  with  the  N+  density  at  30  km  altitude.  These  results 

were  for  the  PCA  Aug  72  day  time  test  profiles  1-3  (see  section  5.1). 

When  data  for  the  PCA  Nov  60  and  PCA  Nov  69,  with  energy  spectra  differing 

from  the  August  72  PCA,  are  included,  this  relation  becomes  more  diffuse, 

cf  Figure  5.9.  Taking  a  variety  of  PCAs  having  differing  spectral  hardness, 

one  can  expect  to  observe  or  predict  a  range  of  field  strengths  for  a 

given  ion  density  at  30  km.  This  is  qualitatively  indicated  in  the  figure 

by  the  shaded  area.  The  unusually  large  signal  attenuation  for  the  PCA  event 

of  12  November  1960  is  due  to  the  unusually  large  ion  pair  production  below 

30  km,  shown  in  Figure  5.3,  which  led  to  enhanced  effects  by  the  resulting 
ions . 
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Figure  5.9  Absolute  field  strength  for  the  WTF  to  Tromso  path  versus  ion 
density  at  30  km  altitude  for  three  PCA  daytime  test  cases 
(Section  5.1),  for  the  PCA’s  of  Nov  60  and  Nov  69,  and  for 
standard  ambient  daytime. 
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6.0  COMPARISON  OF  PCA  SIMULATION  RESULTS 

The  PCA  profiles  shown  in  Figures  5.1  and  5.2  were  used  in  an  inter¬ 
comparison  of  ELF  calculations  over  the  WTF  to  Tromso  path  by  three  laboratories: 
LMSC,  Naval  Ocean  Systems  Center  (NOSC)  and  Pacific  Sierra  Research  (PSR) .  The 
results  show  a  measure  of  agreement  in  some  cases  but  more  often  significant 
disagreements;  see  Table  6.1 

Table  6.1  Computed  change  in  ELF  field  strength  in  going  from  ambient  day 

and  ambient  night  conditions  to  PCA  conditions.  The  calculations 
refer  to  transmission  from  WTF  to  Tromso  at  75  Hz. 


LMSC 

NOSC 

PSR 

PCA  Day  It 1 

-5.7  dB 

-3.8  dB 

-5.6  dB 

It  2 

-2.8 

-1.0 

-3.9 

It  3 

-2.3 

-1.3 

-3.7 

PCA  Night  It  1 

-6.4 

-6.8 

— 

112 

-7.9 

-6.2 

— 

The  differences  may  be  due  to  several  factors: 

In  running  a  code  with  elaborate  segmentation  of  the  path  it  is 
difficult  to  avoid  introducing  differences,  e.g.: 
the  actual  number  of  segments  used  may  differ 
the  positions  and  lengths  of  the  individual  segments  may  vary 

values  of  the  ground  conductivities  may  be  chosen  differently 
because  of  lack  of  agreement  as  to  their  correct  values, 
the  size  assumed  for  the  WTF  area  may  not  be  exactly  the  same 
the  Ne,  N  profiles  supplied  by  LMSC  needed  some  extrapolation  to 
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F  region  altitudes  (PCA  night  #2)  and  to  very  low  altitudes. 

These  extrapolations  may  introduce  minor  differences. 

It  appears  probably  that  a  major  item  involves  the  different  sizes 

of  the  WTF  region  used  by  LMSC  and  NOSC,  which  are  200  km  and  500  km, 

respectively.  The  WTF  transmitter  region  has  low  ground  conductivity 
-4 

at  75  Hz:  3.2  x  10  mho/m.  The  extent  of  this  region  is  taken  to  be  200  km 

in  our  segmentation  of  the  ELF  paths.  W.  Moler,  NOSC,  used  500  km  for  the 

transmitter  region.  However,  at  a  distance  of  500  km  from  WTF  part  of  Lake 
Superior  (a  %  4  ohn/m)  should  be  included.  Thus  the  analyst  faces  a  dilemma 
in  choosing  the  size  of  the  first  segment. 

One  run  was  made  for  PCA  night  112  for  which  we  extended  the  first 
segment  out  to  900  km.  The  results  were: 

PCA  night  112  (Segment  1  size  200  km)  -  162.6  dB 

PCA  night  If 2  (Segment  1  size  900  km)  -  161.5  dB 

It  is  thus  obvious  that  the  size  of  the  segments  for  the  WTF  area  is 
of  importance,  especially  when  one  tries  to  compare  results  from  various  codes 
in  absolute  terms. 
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7.0  EFFECTS  OF  CHANGE  IN  SOLAR  ZENITH  ANGLE  ALONG  THE  PROPAGATION  PATH 

Over  a  long  propagation  path  the  solar  zenith  angle  (x)  will  change 
according  to  the  local  time,  as  will  the  ionospheric  conditions  since  the 
effective  recombination  rates  for  free  electrons  are  dependent  upon  time 
of  day.  Thus,  even  during  a  PCA  when  the  ionization  can  be  assumed  to  be 
constant  over  the  whole  disturbed  polar  cap  region,  the  electron  density 
height  profile  will  vary  with  position. 

In  order  to  estimate  the  magnitude  of  this  effect  during  a  PCA  the 

4  Aug  72  1144  UT-  ionization  profile  was  used  to  calculate  the  N^  profiles 

at  different  solar  zenith  angles,  of  Figure  7.1.  No  change  was  assumed  to 

take  place  in  the  N+  profile.  Waveguide  calculations  were  made  for  several 

simulated  "seasons",  as  exemplified  in  Table  7.1  for  a  spring  equinox  situation 

(21  March  2200  UT)  with  propagation  from  day  conditions  over  the  transmitter 

(X  =  68.5°)  into  night  with  x  =  110°  at  the  receiver  (Tromso,  Norway). 

Runs  were  also  made  with  either  day  or  night  conditions  over  the  complete 

path.  The  results  are  shown  in  Table  7.2,  and  indicate  that  for  more 

accurate  estimates  of  signal  strength  for  a  path  like  WTF  to  Tromso  the 

effect  of  the  solar  zenith  angle  variation  should  be  included. 

Table  7.1  Approximate  conditions  along  the  transmission  path  from  WTF  to 
Tromso  on  21  March  at  2200  UT. 


Segment 

Range  from 
WTF  (km) 

Latitude 

North 

Longitude 

East 

Solar  Zenith 
Angle 

1 

0  -  200 

48° 

271° 

68.  5° 

2 

200  -  1000 

52° 

274° 

72° 

3 

1000  -  1800 

57° 

279° 

77° 

4 

1800  -  2900 

64° 

290° 

8  3° 

5 

2900  -  3300 

68° 

302° 

89.5° 

6 

3300  -  4400 

71° 

318° 

94° 

7 

4400  -  5800 

72° 

355° 

104° 

8 

5800  -  6000 

69° 

18° 

110° 
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ire  7.1  Electron  density  profiles  for  the  4  Aug  72  (1144  UT)  sol 
particle  event  at  the  four  indicated  solar  zenith  angles 
positive  ion  density,  taken  to  be  independent  of  solar  z 
ancle,  is  also  shown. 


Solar 

Zenith 

Path 

Segment 

79.5° 

Whole  path 

110.° 

Whole  path 

79.  5° 

0  -  2900  km 

QO 

CD 

cn 

o 

2900  -  3300 

95° 

3300  -  4400 

110° 

4400  -  6000 

Calculated  Field  Strength 
(75  Hz  (dB  wrt  lA/m) 


156.6 

159.6 


-158.1 


Calculated  Change  from 
Ambient  Night 
(dB  wrt  1  A/m) 


-  1.9 

-  4.9 


-  3.4 
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8.0  ELF  PROPAGATION  CHARACTERISTICS  -  PHYSICAL  INTERPRETATION 

As  described  in  Section  2  important  recent  works  have  shed  light  on 
the  physical  understanding  of  the  ELF  propagation  (Pappert,  1978;  Greifinger 
and  Greifinger,  1978,  1979;  Booker,  1980).  In  this  section  we  shall  describe 
the  Greifinger  and  Greifinger  (1978  and  1979)  work  which  lends  itself 
easily  to  improved  understanding  of  the  ionospheric  parameters  which  are 
important  for  ELF  propagation. 

8.1  The  Greifinger  and  Greifinger  Approach 

Figures  8.1  and  8.2  (Greifinger  and  Greifinger,  1978)  show  the  relation¬ 
ship  between  the  two  heights  h^  and  h^  (see  Section  2),  and  the  variation 
of  the  ELF  wave  fields  and  the  Joule  heating  height  profiles,  respectively. 
The  two  corresponding  scale  heights  G^  and  G^  refer  to  the  scale  heights  at 
these  heights  of  the  conductivity  profile  and  refractive  index,  respectively. 


Figure  8.1  Field  strength  as  a  function  of  altitude  at  a  frequency  of  50  Hz 

for  an  exponential  profile  with  h  =  45  km,  G  =  3  km.  The  fields 
are  normalized  to  |cB$|  -  1  at  the  ground  (Greifinger  and 
Greifinger,  1978). 
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dissipation  (Greifinger  and  Griefinger,  1978). 


of 

km. 

maximum 


It  is  noticed  that  the  ohmic  heating  maximizes  at  two  heights:  h 
ana  slightly  beiow  h.^.  The  ELF  wavefields  rapidly  become  very  small  above 
h^ ,  which  may  then  be  regarded  as  a  "practical"  upper  boundary  of  the  iono¬ 
sphere  as  far  as  ELF  propagation  is  concerned. 

Greifinger  and  Greifinger  (1978)  give  on  the  basis  of  their  theory  the 
following  simple  expressions  for  the  ELF  propagation  parameters: 
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where  S  is  the  eigenvalue  parameter,  v/c  is  the  relative  wave  phase  velocity 
a  is  the  wave  attenuation  coefficient.  This  theory  applies  to  cases  where 
there  is  no  significant  penetration  of  the  electromagnetic  field  to  altitudes 
where  anisotropy  due  to  the  earth's  magnetic  field  must  be  taken  into 
account . 

In  their  1979  paper  Greifinger  and  Greifinger  extend  the  1978  work  to 
include  the  effects  of  anisotropy,  assuming  for  simplicity  a  vertical  terres¬ 
trial  magnetic  field.  The  validity  is  thus  limited  to  higher  magnetic 
latitudes.  The  parameters  which  determine  the  propagation  constant  differ 
for  daytime  and  nighttime  ionospheric  conditions.  In  both  cases,  two 
parameters  which  enter  are  the  frequency-dependent  altitude  hQ,  at  which 
the  conduction  current  parallel  to  the  magnetic  field  becomes  equal  to  the 
displacement  current,  and  the  local  scale  height  of  the  parallel  conductivity, 

a  . 

o 

Under  daytime  conditions,  two  additional  pairs  of  altitudes  (frequency 
dependent)  and  scale  heights  appear  as  parameters.  One  pair  is  the 
altitude  at  which  the  local  reciprocal  wave  number  for  vertically  propagating 
0  waves  becomes  equal  to  the  local  scale  height  of  the  refractive  index, 
and  the  scale  height  of  this  refractive  index.  The  other  pair  is  the 
corresponding  quantities  for  vertically  propagating  X  (whistler)  waves. 

It  is  assumed  in  this  paper  that  these  altitudes  are  attained  in  a  region 
of  the  ionosphere  where  |o  |  >  |a_|,  ou  being  the  Hall  conductivity  and  a 
the  Pedersen  conductivity.  Under  these  conditions,  which  apply  over  a 
substantial  altitude  range,  the  two  pairs  of  parameters  become  identical, 
and  the  analysis  is  somewhat  simpler.  The  vicinity  of  the  single  altitude  h^ 
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is  in  this  case  a  region  of  reflection  rather  than  of  significant  heating, 
as  was  the  case  for  the  isotropic  ionosphere.  The  0  wave  undergoes  nearly 
total  reflection  at  this  altitude,  but  some  of  the  X  wave  energy  leaks 
out  of  the  waveguide,  thereby  contributing  to  the  attenuation. 

For  typical  nighttime  conditions,  a  sharp  reflecting  E-region  bottom 
may  be  encountered  before  the  local  reciprocal  wave  number  becomes  equal 
to  the  local  scale  height.  Under  such  circumstances  the  altitude  of  the 
E-region  bottom  replaces  h^as  a  parameter,  and  the  refractive  index  scale 
height  is  replaced  by  the  local  wavelength  on  the  E-region  side  of  the 
bottom.  Approximate  expressions  for  the  propagation  parameters  are 
derived  similar  in  form  to  those  given  above  for  the  19/8  theory. 

These  theories  have  been  coded  and  have  been  satisfactorily  tested 
at  50  Hz  using  the  exponential  conductivity  profiles  for  which  results  are 
given  in  the  two  papers.  In  addition  calculations  made  at  LMSC  with  the 
NOSC  full  wave  code  using  the  exponential  conductivity  profile  specified 
in  the  1979  paper  yield  results  in  agreement  with  those  obtained  by 
William  Moler  of  NOSC,  which  are  given  in  that  paper. 

In  Figures  8.3  through  8.7,  some  results  of  calculations  applying  the 
Greifinger  and  Greifinger  theory  to  five  of  our  test  profiles  are  given. 

Altitude  profiles  are  shown  of  electron  and  ion  densities,  and 
electron  and  ion  conductivities  and  wavelength  in  the  ionized  medium  at  75  Hz. 

In  Section  9  we  shall  give  some  comparison  between  fullwave  results 
using  the  NOSC  code  and  values  based  upon  the  theory  described  here. 
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8 . 2  Absorpt  ion  and  Re f  lec ti  ion  of  ELK  Waves 

Using  tilt*  PCA  day  and  nighttime  test  cases  supplied  by  LMSC,  W.  Moler 
at  NO SC  has  calculated  Joule  heating  as  well  as  the  altitude  variations  of 
the  integrated  absorption  and  the  reflection  density  of  the  downgoing  wave 
for  these  conditions.  When  compared  with  the  actual  electron  and  ion  density 
variation  with  height,  the  variations  of  the  above  entities  reveal  interesting 
characteristics  of  ELF  propagation. 

As  an  example,  in  Figures  8.8  and  8.9  we  have  plotted  the  profiles 
for  the  PCA  day  test  case  //I  for  75  Hz.  Figure  8.8  clearly  shows  the 
importance  of  the  enhanced  ion  concentration,  peaking  approximately  at  40  km. 
The  relative  heating  loss  (dotted  curve)  indicates  that  the  maximum  Joule 
heating  occurs  at  35  km.  Integrated  absorption  in  dB/1000  km  (dash-dotted 
curve)  furthermore  shows  that  approximately  1.3  dB  out  of  the  total  2.1 
dB/1000  km  can  be  associated  with  this  layer.  The  free  electrons  cause 
Joule  heating  mainly  at  higher  altitudes,  maximizing  around  60  km.  In 
Figure  8.4  the  reflection  density  profile  of  the  downgoing  wave  (dash- 
dotted  curve)  indicates  that  the  region  between  50  and  80  km  is  the  main 
reflecting  region  during  the  PCA  conditions. 

Similar  curves  have  been  obtained  for  other  and  N+  profiles.  In 
the  case  of  quiet  ambient  nighttime  conditions,  the  ELF  fields  penetrate 
to  greater  heights  before  reflection  takes  place.  In  this  case  the  Joule 
heating,  essentially  due  to  electrons,  occurs  within  quite  a  narrow  altitude 
interval  between  65  and  75  km,  whereas  the  ions  cause  another  maximum  at 
around  110  km.  Also,  in  this  case  there  is  only  an  insignificant  contribution 
to  the  attenuation  from  Joule  heating  below  40  km.  Significant  contributions 
to  tiie  H,  field  (at  ground  level)  are  coming  from  altitudes  between  100-120  km 
but  altitudes  up  to  about  160  km  do  contribute. 
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ABSORPTION  (dB/1000  km)  — 


PC  A  DAY  CASE  NO.  1 


ELECTRON,  ION  DENSITY  (cm  ) 


RELATIVE  HEATING  LOSS/UNIT  VOLUME . 

Figure  8.8  Altitude  variation  of  ionospheric  and  ELF  characteristics 
for  PCA  daytime  conditions  (test  case  #1). 

-  N+,  ion  density  (heavy,  solid  curve) 

-  Ne,  electron  density  (heavy,  dashed  curve) 

-  Relative  heating  loss/unit  volume  (dotted  curve) 

-  Integrated  Absorption  (dB/1000  km)  (dash-dotted  curve) 

The  last  two  curves  were  computed  by  W.  Moler,  NOSC,  for  a  frequency 
of  75  Hz  for  the  pictures  N+  and  Ng  profiles  (Moler,  1978). 
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9.0  ACCURACY  ASSESSMENT  OF  ELF  CALCULATIONS 


In  discussing  the  accuracy  of  the  ELF  calculations  one  should  keep 
in  mind  the  distinction  between  accuracies  (or  inaccuracies)  inherent  in 
the  actual  mathematical  treatment  of  the  propagation  problem  and  effects 
due  to  less  than  complete  knowledge  of  the  propagation  medium  and  its 
spatial  and  temporal  irregularities,  the  ionospheric  conductivity  profile 
etc.  In  an  actual  ELF  calculation,  however,  it  may  not  be  possible  to  make 
a  sharp  distinction  between  these  two  classes  of  possible  errors. 

As  as  example  we  show  in  Table  9.1  an  intercomparison  of  predicted 

ELF  field  strengths  for  ambient  and  disturbed  condtions  based  upon  use  of 

the  identical  mathematical  computer  code  (the  NOSC  code).  The  discrepancies 

(in  most  cases  less  than  1  dB)  are  attributed  to  slightly  different  ways 

of  segmenting  the  propagation  path,  different  values  for  ground  conductivities, 

or  other  small  differences  in  organizing  the  computer  code  run. 

Table  9.1  Comparison  between  ELF  absolute  field  strength  calculations  at 
75  Hz  for  the  path  WTF  to  Tromso. 


Condition 

LMSC 

dB/lA/m 

NOSC 

dB/lA/m 

Day,  ambient 
Night,  ambient 

-154.7 

-154.7 

-155.4 

-154.6 

Day  PCA  //I 

Day  PCA  it  2 

Day  PCA  it  3 

-160.4 

-157.5 

-157.0 

-159.2 

-156.4 

-156.7 

Night  PCA  //I 
Night  PCA  it 2 

-161.1 

-162.6 

-161.4 

-160.8 

*  The  NOSC  data  supplied  by  W.  Moler  have  been  adjusted  by  -1.6dB  due  to  a 
different  value  used  for  the  soil  conductivity  under  the  WTF  transmitter 
(NOSC  data  were  actually  computed  for  a  =  2.2  10“^  mho/m,  whereas  the 
field  strength  values  in  the  table  are  normalized  to  o  =  3.2  10~^  mho/m, 
which  is  used  by  LMSC). 
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Different  mathematical  approaches  to  the  propagation  problem 
may  yield  results  with  less  than  perfect  agreement.  As  an  example,  the 
results  of  a  comparison  between  altitude  profiles  of  the  Joule  heating 
calculated  by  K.  C.  Field  and  by  W.  Moler  for  PCA  day  test  case  i/1  is 
given  in  Figure  9.1.  The  curves  have  been  normalized  to  provide  equal 
values  at  the  lower  maximum  in  the  30-35  km  altitude  regime  and  they  indicate 
that  a  reasonable  agreement  exists  between  two  approaches.  Both  calculations 
show  that  tlie  largest  attenuation  of  75  Hz  waves  occurs  between  approximately 
25  to  AO  km  with  a  secondary  maximum  around  60  km.  The  former  maximum 
may  be  linked  with  high  ion  densities,  whereas  the  latter  maximum  is  mainly 
due  to  the  free  electrons. 

Such  differences  may  indicate  the  uncertainty  inherent  in  the  various 
theoretical  approaches. 

A  further  example  of  this  point  is  provided  in  Table  9.2,  where  we 
compare  results  obtained  for  a  variety  of  electron/ion  profiles  by  use 
of  the  full  wave  method  and  the  Greifinger  and  Greifinger  approach  discussed 
in  Sections  2  and  8.  It  is  noticed  that  some  differences  are  present, 
but  the  remarkable  agreement  between  the  full  wave  method  and  a  simplified 
theory  is  striking.  In  consideration  of  Table  9.2,  it  should  be  noted 
that  under  highly  disturbed  conditions  the  approach  used  in  the  Greifinger 
and  Greifinger  (1978)  theory  is  preferred  as  stated  in  their  1979  paper. 
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RELATIVE  HEATING  LOSS  /  UNIT  VOLUME  (ARB.  SCALE) 


Figure  9.1  Comparison  of  Joule  heating  for  PCA  day  test  case  #1 
by  E.  C.  Field  (PSR)  and  W.  Moler  (NOSC) .  The  curves 
have  been  normalized  to  give  equal  values  at  the  lower  maxi 
mum. 
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Table  9.2  Comparison  oi  calculation  of  KLF  relative  phase  velocity  and 

attenuation  rate  af  75  !iz  using  the  fall  wave  method  (LMSC)  and  the 
approximate  theory  of  Greifinger  &  Greifinger  (1978,  1979).  LtlSC 
values  are  given  for  the  WTF  area  and  the  Norwegian  Sea. 


Ambient 

Night 


G&G  78 
79 
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10.0  AMBIENT  PROFILES  OF  ELECTRON  AND  ION  DENSITY 
10.1  Introduction 

Although  our  studies  of  ELF  propagation  are  mostly  concerned  with  the 
ionosphere  perturbed  by  precipitating  particles  it  has  been  necessary  to  use 
ambient  profiles  of  electron  and  ion  density.  Our  long  path  calculations, 
e.g.  from  WTF  to  Tromso,  are  segmented  following  ground  conductivity  varia¬ 
tions  and  frequently  one  or  more  segments,  such  as  the  first  one  at  the  WTF 
transmitter  are  considered  to  have  ambient  ionospheric  conditions.  For 
comparison  purposes,  calculations  have  also  been  performed  in  which  all  seg¬ 
ments  of  the  long  parth  are  given  the  same  ambient  profile.  The  profiles 
were  typically  those  used  by  Pappert  and  Moler  and  colleagues  (e.g.  Pappert 
and  Moler,  1974).  These  profiles,  shown  in  Fig.  3.2,  have  come  to  be  called 
the  "standard"  ambient  day  and  night  profiles. 

For  modeling  the  ionosphere  realistically  over  a  long  segmented 
path,  taking  into  account  the  variations  along  the  path  in  local  time  as 
well  as  atmospheric  density  and  temperature,  a  good  ion  chemistry  model  is 
needed.  Ideally,  the  model  should  compute  accurately  ambient  as  well  as 
disturbed  conditions.  Our  model  has  been  adjusted  to  fit  PCA  daytime 
disturbed  conditions  during  the  2  November  1969  and  4  August  1972 
events.  The  model  has  been  used  in  an  effort  to  make  realistic  calculations 
of  those  events  including  ambient  conditions  over  the  WTF  transmitter. 

Ambient  calculations  foi  >  1972  event  are  shown  in  Figure  10.1.  The 
figure  compares  the  standard  ambient  day  profiles  with  those  calculated  for 
WTF  at  one  time  during  the  actual  event  and  during  a  simulated  event  such  as 
the  4  August  1972  event  but  shifted  in  season  to  the  winter  solstice.  It  may 
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Figure  10.1  Comparison  of  ambient  daytime  electron  density  profiles  at  WTF 
for  the  two  seasons  and  times  of  day  with  the  standard  ambient 
daytime  profile  discussed  earlier  in  Section  3. 
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be  noted  that  these  calculated  profiles  differ  from  the  standard  profile, 
having  larger  electron  densities  at  the  lower  altitudes.  As  will  be  shown 
later,  recent  daytime  measurements  of  ambient  profiles  tend  to  confirm  such 
larger  densities.  These  larger  densities  are  significant,  leading  to  a 
weakening  of  the  ELF  signal  at  75hz  over  the  WTF  to  Tromso  path  by  0.7  dB 
relative  to  standard  ambient  day  used  over  the  WTF  segment. 

There  have  been  many  measurements  of  electron  and  ion  density  profiles 
since  the  "standard"  profiles  of  Figure  3.2  were  generated,  particularly 
during  daytime.  There  are  fewer  nighttime  measurements  and  because  of  their 
scatter  it  would  be  difficult  to  improve  the  present  standard  night 
profiles.  Also  the  present  nighttime  profile,  used  in  NOSC  code  calculations 
for  the  WTF  to  Connecticut  path  yielded  good  agreement  with  measurements 
made  over  the  path  (Imhof  et  al.,  1976).  The  purpose  of  this  section  is  to 
propose  a  new  standard  daytime  profile  based  on  recent  measurements. 

10.2  Proposed  New  Standard  Ambient  Profile 

In  developing  this  profile  considerable  use  was  made  of  the  work  of 
Keefer  (1975)  who  assembled  some  representative  measured  ambient  profiles. 
Several  daytime  profiles  are  shown  in  Figure  10.2  along  with  the  present 
standard  daytime  profile.  All  of  the  experimental  profiles  are  based  on 
mid-latitude  measurements  made  at  a  solar  zenith  angle  near  x  =  60°.  Each 
of  the  three  profiles  within  the  range  60  to  160  km  represents  a  concensus 
by  the  authors  of  a  number  of  measurements.  The  E  region  profile  of  Maeda 
(1972)  was  based  on  winter,  low  solar  activity  conditions  and  differs  very 
little  from  the  standard  ambient  day  model.  The  Mechtly  et  al.  (1972)  curve 
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is  a  mean  of  quiet  sun  measurements.  The  curve  of  Mitra  and  Chakrabarty 
(1971)  is  a  concensus  of  measurements  made  by  a  number  of  different  methods 
in  the  65-90  km  range.  The  lower  altitude  profile  is  based  on  the  work  of 
York  (e.g.,  1980)  who  has  been  reevaluating  rocket-borne  probing  techniques 
for  measuring  electron  densities,  among  them  the  blunt  probe  (York,  1979). 

This  re-analysis  leads  to  considerably  larger  electron  densities  below  60  km 
than  specified  by  the  present  standard  ambient  day  profile.  We  have  based 
the  profile  shown  between  40  and  70  km  on  measurements  reported  by  York  (1980) 
made  at  White  Sands  in  October  at  a  solar  zenith  angle  near  60°. 

These  more  recent  results,  as  shown  in  Figure  10.2,  suggest  that  the 
present  daytime  standard  ambient  profile  is  considerably  low  at  altitudes 
below  60  km,  is  somewhat  too  high  between  60  and  100  km,  and  above  100  km  is 
reasonably  good.  We  suggest  that  the  new  composite  ambient  profile  shown  in 
Figure  10.3  be  used  for  comparison  purposes  to  represent  daytime  conditions 
in  future  ELF  propagation  studies.  Of  course,  ambient  electron  density  profiles 
depend  on  solar  zenith  angle,  latitude,  season  and  solar  activity  and  no  single 
profile  can  represent  all  of  these  conditons.  Also  shown  in  Figure  10.3 
are  the  present  standard  profile  and  a  profile  (x  =  75°)  used  by  Booker  (1980). 

Figure  10.3  also  suggests  a  positive  ion  ambient  profile  to  be  used  with 
the  new  electron  profile.  It  was  calculated  using  relations  from  simple  lumped 
parameter  ionization  theory: 

N+  =  Ne  +  N~,  (1) 

specifying  charge  neutrality,  and 

Q  =  cl  N+N  +  aTN+N_  (2) 

^  D  e  X 
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the  steady  state  equation  for  positive  ions,  where  Q  is  the  ion  pair 
production  rate,  is  the  electron-ion  recombination  rate,  and  ot^.  is  the 

ion  neutralization  rate.  Equations  (1)  and  (2)  may  be  combined  to  yield 

N+  =  (/  x2  + y  -  x)/2  (3) 

_6  2  —X 

where  x  =  -  1  )  and  y  =  4  Q/a^.  The  value  3  x  10  cm  sec 

was  assumed  for  and  the  values  of  are  those  of  Table  5.1  The  values  of 
Q,  galactic  cosmic  ray  (GCR)  production  rates,  are  calculated  from  the  formula 
of  Heaps  (1978)  for  average  solar  conditions  and  a  magnetic  latitude  of  57°  N 
corresponding  approximately  to  the  Wisconsin  Test  Facility  location.  The 
GCR  Q  rates  were  augmented  with  Lyman  alpha  production  rates  above  70  km. 

Calculations  have  been  made  of  ELF  propagation  at  75  Hz  from  WTF  to 
Tromso  using  the  suggested  new  profiles  supplemented  by  the  old  standard  values 
below  40  km  for  comparison  with  similar  calculations  for  the  old  standard. 

Of  course  the  segmentation  and  ground  conductivities  were  identical.  The 
resulting  calculated  overall  signal  attenuation  is  1.0  dB  larger  over  the 
6000  km  path,  a  significant  change.  This  change  is  due  to  the  larger  electron 
densities  in  the  new  profile  near  50  km.  The  effect  of  the  new  profile  is  to 
produce  larger  excitation  factors  which  tend  to  increase  the  received  signal 
but  larger  attenuations  which  reduce  it  to  produce  the  net  effect  at  the 
receiver  of  signal  strength  1  dB  smaller. 
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11.0  ION  CHEMISTRY  MODEL  IMPROVEMENT 
11.1  Introduction 

Because  of  the  strong  dependence  of  ELF  propagation  on  ionospheric  con¬ 
ditions  demonstrated  in  this  and  other  studies  is  it  very  important  to  have 
the  capability  of  making  accurate  predictions  of  electron  and  ion  density 
profiles,  given  altitude  profiles  of  ion  pair  production  rates  derived  from 
satellite  or  rocket  measurements  of  energetic  particle  precipitation. 

This  is  particularly  important  for  disturbed  conditions  such  as  those 
occurring  in  polar  cap  absorption  events  because  of  the  potentially 
large  effects  in  reducing  ELF  signals  transmitted  through  the  polar  regions. 
The  theory  of  ELF  propagation  in  the  earth-ionosphere  waveguide  has  been 
advanced  as  discussed  earlier  in  this  report  to  the  point  where  Booker  (1980) 
states  that  "at  present,  and  probably  for  a  long  time  to  come,  overall  error 
in  such  calculations*  is  likely  to  be  controlled  by  imperfect  knowledge  con¬ 
cerning  worldwide  behavior  of  the  ionosphere".  The  pertinent  available  data 
on  the  ionosphere  should  be  used  to  improve  the  predictive  capability  of 
ionospheric  models.  A  good  model  will  be  needed  when  a  field  campaign  is 
conducted  to  test  ELF  propagation  theory  during  a  PCA  event.  In  order  to 
validate  ELF  propagation  theory,  over  long  test  paths,  ionospheric  prediction 
will  be  needed  to  extend  the  ionospheric  measurements  which  can,  as  a  practice 
matter,  be  made  only  over  limited  time  spans  and  limited  spatial  extents. 

One  of  our  present  objectives  is  to  use  more  of  the  electron  density  data 
measured  during  the  4  August  1972  solar  particle  or  PCA  event  by  the  SRI 
Chatanika  radar  than  has  been  done  in  our  previous  modeling  efforts.  As 
discussed  in  our  previous  report  (Imhof  et  al.  1980)  previous  modeling  has 


♦Calculations  such  as  those  made  by  Booker  (1980)  and  by  use  of  the  NOSC  code 
(Pappert  and  Moler  1974,  1978). 
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been  based  on  data  at  some  28  points  in  time.  The  ion  chemistry  model  des¬ 
cribed  in  that  earlier  report  agrees  reasonably  well  with  the  daytime  data 
from  the  August  1972  event,  and  also  with  measurements  made  during  the 
2  November  1969  PCA  (Swider  et  al.  1978).  These  results  are  shown  in  Figures 

11.1  and  11.2.  There  was  however,  a  lack  of  agreement  in  magnitude  with  the 
measurement  made  at  one  time  at  60  km  altitude  in  deep  twilight  near  midnight 
August  4/5  1972;  the  model  predicted  a  significantly  smaller  value.  The  data 
from  the  1969  event  measured  during  sunrise  show  a  time  lag  in  the  rise  of 
electron  density  relative  to  that  predicted  by  our  model  and  by  the  model  of 
Swider  et  al.  (1978).  We  have  devoted  our  attention  particularly  to  these 
twilight  modeling  problems  during  the  present  effort.  Accurate  modeling  of 
twilight  effects  is  important  for  prediction  of  ELF  propagation  through  twi¬ 
light  regions  near  the  solar  terminator. 

11.2  The  SRI  Data 

These  data,  made  available  by  R.  Vondrak  of  SRI  International,  consist  of 
raw  electron  densities  measured  by  the  Chatanika  incoherent  scatter  radar  in 
the  time  period  3-7  August  1972  UT  including  the  main  peak  of  the  event  on 

4  August.  The  data  are  based  on  10  minute  averages  of  the  received  signal  and 
are  available  typically  for  over  100  of  the  144  such  periods  during  the  day. 

The  data  has  been  interpolated  between  the  actual  altitudes  of  measurement  to 
altitudes  from  50  to  80  km  separated  by  5  km.  The  data  provided  was  not  fully 
corrected  for  receiver  recovery  and  a  second  correction  has  been  applied 
following  the  analysis  of  Watt  (1975).  The  data  for  70,  60  and  50  km  on  4  and 

5  August  1972  UT  are  presented  in  Figures  11.3  to  11.5.  The  60  and  70  km  data  are 
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is  /Q  multiplied  by  a  normalizing  factor.  The  solid  line  represents  electron 
density  calculated  by  our  ion  chemistry  model.  See  test  for  further  explanation. 
The  receiver  recovery  correction  deduced  by  Watt  (1975)  is  shown. 
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density  calculated  by  our  ion  chemistry  model.  See  test  for  turther  explanation 
The  receiver  recovery  correction  deduced  by  Watt  (1975)  is  shown. 
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10  minute  averages  and  the  SO  km  data,  because  of  the  large  scatter  have  been 
converted  essentially  to  30  minute  averages  by  averaging  each  original  point 
where  possible  with  its  neighbors.  The  experimental  electron  density  data 
have  not  been  corrected  for  the  presence  of  negative  ions.  The  effect  of 
negative  ions  is  to  require  a  reduction  of  the  electron  densities  by  a  factor 
ranging  between  one  and  two,  depending  on  the  negative  ion-electron  density 
ratio  as  discussed  in  our  previous  report  (Imhof  et  al.,  1980).  The  result  of 
this  correction  is  to  reduce  nearly  all  of  the  data  at  50  and  60  km  by  a 
factor  of  two  and  to  significantly  reduce  the  values  at  70  km.  There  is  also 
plotted  on  each  figure  values  of  the  square  root  of  the  ion  pair  production 
rates,  /Q,  derived  from  measurements  of  precipitating  proton  fluxes  and 
energies  made  by  Lockheed  instrumentation  on  polar-orbiting  satellite 
1971-089A  (Reagan  et  al.,  1980).  The  values  of  /Q  have  been  multiplied  by  a 
normalizing  factor,  different  for  each  altitude  so  as  to  make  them  comparable 
to  the  electron  densities. 

The  rationale  behind  such  use  of  the  /Q  is  based  on  the  expression  some¬ 
times  used  in  ionospheric  analyses: 

Q  “  a  „  N  2 
ef  f  e 

where  is  an  effective  loss  coefficient  for  electrons.  Then  the  electron 

density  may  be  expressed  as 

s«  ■ 
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so  that  if  aeff  were  constant  the  time  dependences  of  Ng  and  /Q  would  be  the 
same.  Expressed  in  terms  of  lumped  parameters 

«eff  “  (l+XXop+Xoj  ) 

where  X  is  N~/Ne,  the  negative  ion  to  electron  ratio,  is  the  effective 
electron-positive  ion  dissociative  recombination  coefficient,  and  oij  is  the 
effective  positive  ion-negative  neutralization  coefficient.  As  implied  a 
lumping  together  of  all  the  positive  ions  and  of  all  the  negative  ions  is 
implied  by  this  simplified  expression.  Although  aj  is  thought  to  be  rela¬ 
tively  constant,  aD  and  particularly  X  vary  with  minor  species  concentrations 
and  with  photoeffects  and  are  thus  time  dependent.  Thus  Ng  is  not  expected  to 
show  the  same  time  dependence  as  /Q.  However,  a  comparison  of  the  two  does 
show,  for  example  that  the  main  peak  of  the  event  probably  is  larger  than  the 
Q  values  would  indicate.  The  Q  values  are  grouped  in  pairs  separated  by  about 
1.6  hours,  the  satellite  orbital  period. 

The  electron  densities  computed  by  our  ion  chemistry  model  (Imhof  et  al., 
1980)  are  also  plotted  on  the  figures.  That  model,  computes  the  diurnal 
variation  of  minor  neutral  species  such  as  0,  A),  O3,  NO,  and  NO2  which 

have  important  influences  on  the  ion  chemistry.  The  ion  chemistry  model  which 
is  coupled  to  the  neutral  chemistry  is  a  simplified  model.  In  this  simplifi¬ 
cation,  for  example,  all  of  the  hydrated  protons  H+(HjO)n  (n>0)  are  lumped 
into  one  species,  as  are  the  hydrates  of  N0+.  There  are  four  negative  ions: 
02~,  X~,  Y~,  Z~.  X-  represents  a  lumping  of  0^~  and  C0^_;  Y~  represents  0~, 
O3-  and  CO^-;  and  Z~  represents  O2TNO,  NC^-  and  NO3-.  The  negative  ion  reac- 
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tion  scheme  is  illustrated  in  Figures  11.6  and  11.7.  The  diurnal  variations 
of  the  photocoefficient  are  computed  in  advance  by  a  separate  code.  The  model 
computes  one  altitude  at  a  time,  i.e.  it  is  a  box  model. 

11.3  The  Twilight  Problems 

In  Figure  11.4  showing  the  SRI  data  for  60  km,  the  time  of  particular 
interest  is  the  twilight  period  around  midnight  4/5  August  local  (Alaska 
Standard)  tims.  At  the  latitude  of  Chatanika,  Alaska  (65.1°N)  at  60  km  alti¬ 
tude  on  4  August  the  sun  does  not  set  completely  at  midnight,  the  solar  zenith 
angle  becomes  no  larger  than  about  98°.  That  altitude  is  still  illuminated 
by  visible  light;  most  of  the  UV  has  been  absorbed  in  the  atmosphere.  It  will 
be  noted  in  the  figure  that  Q  is  relatively  constant  at  this  time.  The  SRI 
data  near  midnight  have  a  mean  value  several  times  larger  than  the  model  (refer 
also  to  Figure  11.1)  and  thus  there  is  a  discrepancy.  It  also  is  evident  that 
the  data  and  the  model  have  different  time  dependences  at  sunrise:  the  data 
appear  to  rise  more  slowly.  The  corresponding  time  dependence  at  the  im¬ 
mediately  preceding  sunset  is  not  clear  partly  because  of  the  rapid  drop  in  Q 
in  the  early  evening  which  may  have  influenced  the  minor  neutral  species  and 
hence  the  ion  chemistry.  There  is  a  question  of  how  well  the  subtraction  of 
the  receiver  recovery  correction  from  the  raw  data  has  been  done.  It  may  be 
noted  at  the  beginning  of  the  60  km  plot  before  the  start  of  the  event, 
experimental  values  of  Ng  tend  toward  zero  on  this  plot  as  would  be  expected 
for  near  ambient  conditions. 

Examination  of  the  50  km  data  in  Figure  11.5  during  the  same  time  period 
around  midnight  4/5  August  local  time  shows  a  larger  discrepancy  between  data 
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Figure  11.6  Schematic  of  part  of  the  known  negative  ion  reaction  scheme 
showing  the  lumping  together  of  ions  used  in  the  simplified 
scheme . 
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and  model.  The  receiver  recovery  correction  made  at  50  km  i6  larger  and  may  be 
a  source  of  some  error.  The  scatter  of  the  50  km  data,  even  with  30  minute 
averaging,  is  too  large  to  permit  a  comparison  of  time  dependences. 

There  is  evidence  about,  the  twilight  time  dependence  of  from  an  ear¬ 
lier  investigation  by  Reid  (1961)  of  enhanced  ionization  produced  ly  solar 
proton  bombardment  during  PCA  events.  Measurements  were  made  by  riometers 
during  several  PCA  events  of  the  vertical  absorption  of  30  and  60  Mhz  cosmic 
noise.  The  absorption  was  large  during  daytime  when  electron  densities  were 
large  but  decreased  at.  night  only  by  a  factor  of  about  A  indicating  qualita¬ 
tively  that  electron  densities  had  decreased  in  approximately  the  same  ratio. 

Reid  used  balloon  measurements  of  the  precipitating  protons  to  deduce  ion 
pair  production  altitude  profiles,  and,  assuming  constant  spatially  uniform 
proton  fluxes  computed  differential  absorption  (db/Km)  altitude  profiles  at  the 
riometer  frequencies.  The  absorption  at  30  Mhz  had  a  fairly  definite  peak  at  an 
altitude  of  about  58  km  although  the  whole  region  from  35  to  80  km  was  effective. 

To  try  to  explain  the  measured  time  dependence  of  the  riometer  absorption 
at  sunset  and  sunrise  and  to  explain  the  residual  nighttime  absorption,  Reid 
set  up  a  simple  ion  chemistry  model  which  included  electron-ion  recombination, 
attachment  to  form  O2  ,  photodetachment  and  collisional  detachment  from  >  and 
ionic  recombination.  Good  data  were  available  for  ionic  recombination  and  for 
the  attachment  to  and  photodetachment  from  07~  which  is  due  largely  to  light  in  the 
visible.  The  estimated  electron-ior.  recombination  rate  was  3  x  10  cmJ 
sec  In  order  to  maintain  the  electron  density  necessary  to  provide  the 

observed  nighttime  absorption  it  was  necessary  to  assume  a  rate  of  collisional 
detachment,  e.g.  O2  +  M  +  0£  +  e  +  M  of  about  10”^  cm^  sec-^.  The  rate  of 
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such  a  reaction  is  known  to  be  several  orders  of  magnitude  smaller,  and  there 
does  not  appear  to  be  according  to  current  knowledge,  any  detachment  process 
which  can  provide  the  explanation  for  the  small  day-night  absorption  ratio. 

The  observed  twilight  time  dependence  of  the  riometer  absorption  also 
did  not  match  the  model.  The  decrease  of  absorption  at  sunset  started  at 
an  earlier  solar  zenith  angle  (x  ^  88°)  than  did  the  model  (x  %  92°)  and 
at  sunrise  the  rise  of  absorption  ended  even  later  (x  ^  85°)  than  the  onset 
of  sunset,  suggesting  an  asymmetry  between  sunset  and  sunrise.  The  suggestion 
was  offered  that  the  time  dependences  might  be  explained  by  detachment  from 
an  unknown  ion  by  UV  radiation. 

Although  we  know  much  more  about  mesospheric  negative  ion  chemistry 
today  than  was  known  in  1961  we  still  cannot  answer  in  detail  the  questions 
raised  by  this  interesting  work  of  Reid. 

Measurements  made  at  night  during  a  PCA  event  on  18  November  1968  over 
Churchill,  Canada  at  0200  CST  have  been  advanced  by  Swider  et  al.  (1971) 
as  evidence  that  nighttime  electron  densities  can  be  explained  by  a  model  in 
which  there  is  no  nighttime  detachment.  They  cite  data  measured  during  a 
nighttime  nuclear  event  (Whitten  et  al.  1965)  as  further  evidence  for  this 
point  of  view.  Also,  some  data  measured  during  the  PCA  of  2  November  1969 
(Ulwick,  1972)  fits  our  model  as  indicated  in  Figure  11.2.  Our  model  has 
no  nighttime  detachment  mechanism  at  that  time  of  year. 

In  summary,  there  is  evidence  from  the  riometer  measurements  analyzed 
by  Reid  (1961)  that  there  may  be  a  detachment  process  effective  in  the  dark 
polar  night  during  a  PCA  event.  There  is  evidence  from  the  SRI  Chatanika 
radar  measurements  made  during  the  August  1972  event  in  deep  twilight  that 
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visible  light  photodetachment  or  some  other  process  may  be  occurring.  On 
the  other  hand,  measurements  made  at  night  during  the  November  1969  PCA 
event  apparently  show  no  detachment.  We  have  been  unable  to  find  any  way 
of  resolving  these  differences. 

There  is  a  lack  of  agreement  between  measured  time  dependences  of 
electron  density  at  twilight  during  PCA  events  and  the  results  of  current 
and  earlier  models.  It  may  however  be  possible  to  shed  some  light  on  this 
question  as  discussed  in  the  next  section. 

11.4  The  Twilight  Time  Dependence 

Swider  et  al.  (1978)  present  sunset  and  sunrise  electron  density 
measurements  (made  by  rocket)  during  the  PCA  event  of  November  1969.  The 
results  clearly  show  a  sunset-sunrise  hysteresis  at  60  and  70  km  altitude 
i.e.  the  rise  of  the  electron  density  at  sunrise  is  not  symmetrical  and 
occurs  later  than  would  be  expected  if  it  were  symmetrical  around  midnight 
with  sunset.  Also  presented  are  results  from  a  time-dependent  ion  chemistry 
model  which  agree  quite  well  with  the  sunset  measurements  but  shows  no 
hysteresis,  i.e.  the  modeled  sunrise  electron  densities  rise  more  rapidly 
than  the  measured  values.  The  model  includes  much  of  the  currently  known 
chemistry  of  negative  ions  in  the  D  region.  It  appears  that  detachment  in 
the  model  is  probably  largely  photodetachment  rather  than  chemical  detachment 
by  0  atoms  or  excited  (^(^A)  molecules. 

Our  box  model  described  above  when  applied  to  either  the  November  1969 
PCA  event  or  the  August  1972  event  during  the  twilight  period  around  midnight 
4/5  August,  also  shows  little  or  no  hysteresis  and  tends  at  60  km  to  be 
dominated  by  photodetachment  in  the  early  stages  of  sunrise  e.g.  at  0300 
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local  time.  Later  in  sunrise  detachment  from  0^  by  0  atoms  and  (  A) 
molecules  becomes  commensurate  with  photodetachment. 

The  time  dependence  given  to  all  of  the  photodetachment  processes 
in  the  box  model,  for  lack  of  further  information,  usually  has  been  that 

O 

computed  for  0  using  the  solar  flux  in  the  3000-8000  A  range.  This  tends 
to  cause  the  electron  density  at  sunrise  to  rise  too  early  and  too  rapidly 
compared  with  measurements.  When  the  photodetachment  in  the  model  from  ions 
other  than  0^  are  given  the  same  time  dependence  as  the  photodissociation 
process . 

03  +  hv  -  0(1D)  +  02(1A)  (1) 

O 

which  occurs  for  wavelengths  shorter  than  about  3100  A,  photodetachment 
starts  later  than  does  the  0  photodetachment  but  rises  too  rapidly. 

For  the  box  model,  the  photocoefficients  and  their  time  dependences 
are  computed  in  advance  using  a  separate  code  having  selected  but  fixed 
altitude  profiles  of  the  absorbing  species  0^  and  0^.  In  reality  the  0^ 
profile  is  time-varying  during  twilight  and  thus  consideration  should  properly 
be  given  to  the  effect  which  the  temporal  dependence  of  the  ozone  may  have 
on  the  photodetachment  coefficients. 

We  have  made  a  preliminary  study  of  this  question  by  making  use  of  a  one¬ 
dimensional  (1-D)  neutral  chemistry  model  developed  several  years  ago  on 
another  project  for  the  purpose  of  finding  the  variations  of  minor  species 
densities  in  the  mesosphere  during  the  August  1972  PCA  event.  It  computes  a 
neutral  chemistry  similar  to  that  of  the  box  model  simultaneously  at  41  equally 
spaced  altitudes  between  40  and  120  km.  The  photocoefficients  are  updated  at 
selected  times  using  the  changing  0^  profile. 
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The  1-D  code  was  first  run  through  several  complete  ambient  days  with 
atmospheric  densities  and  temperatures  approximating  those  of  Chatanika, 

Alaska  (~65°N)  in  August.  Then  it  was  used  to  model  approximately  the 
August  1972  event  up  to  noon  on  5  August.  Ion  pair  production  altitude 
profiles  varying  with  time  were  used  to  augment  the  odd  nitrogen  species 
by  adding  1.5  N  atoms  per  ion  pair,  and  the  odd  hydrogen  production  accom¬ 
panying  H+  (^0)  ion  formation  and  removal  was  simulated  by  adding  one 
H  atom  and  one  HO  radical  for  each  ion  pair  produced,  the  upper  limit.  In 
addition  there  is  some  production  of  0  atoms  accompanying  the  ionization 

The  results  at  60  km  showed  that  despite  the  changing  ozone  used  in  the 
coupled  photocoefficient  code  the  coefficients  for  the  processes  (1)  above, 

03  +  hv  -*  0  +  02  A  >  3100  A°  ,  (2) 

and 

02  +  hv  -►  0  +  0  (3) 

showed  very  little  hysteresis.  The  time  dependences  of  these  processes  were 
slightly  different  from  those  of  the  box  model,  although  the  photocodes  used 
in  the  two  models  are  identical.  The  differences  must  be  due  to  the  differing 
ozone  profiles  input  to  the  photocodes.  The  fact  that  the  photodissociation 
processes  show  little  or  no  hysteresis  probably  indicates  that  photodetachment 
processes  involving  UV  also  do  not. 

The  most  interesting  result  was  that  a  pronounced  hysteresis  appeared  in 
the  0  atom  density  at  60  km  as  shown  in  Figure  11.8.  The  sunrise  lag  in 
the  rise  of  0  atoms  may  be  understood  by  considering  qualitatively  some 
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features  of  the  odd  oxygen  (0,  0^)  chemistry.  From  a  daytime  high  the  0 
atom  density  and  also  the  ozone  density  tend  to  fall  with  the  decrease 
during  sunset  in  the  rate  of  photodissociation  of  C^,  process  (3).  As  the 
rate  of  process  (1)  begins  a  rapid  decrease,  0  atoms  decrease  to  a  low 
level,  reacting  with  0^  to  form  0^  which  increases  to  a  nocturnal  peak. 
Throughout  the  night  the  ozone  is  slowly  removed  by  catalytic  reactions 
involving  odd  nitrogen,  odd  hydrogen  species,  and  a  few  0  atoms.  At 
sunrise,  the  nocturnal  ozone  has  been  depleted  and  because  of  the 
reduced  amount  of  0^  the  production  of  0  atoms  by  processes  (1)  and  (2)  is 
slow  enough  to  cause  the  hysteresis  or  relative  lag  in  the  rise  of  0 
atom  density.  Although  not  specifically  calculated  in  the  1-D  model, 
there  must  be  a  closely  related  hysteresis  in  time  dependence. 

Similar  0  atom  hysteresis  was  also  found  to  occur  for  ambient  twilight 
conditions  a'c  60  km  in  the  1-D  model  and  might  also  be  expected  to  occur  in 
the  box  model.  In  fact,  such  an  effect  does  occur  at  60  km  in  the  box  model, 
but  much  smaller  in  magnitude  than  in  the  1-D  model,  not  large  enough  to 
have  an  appreciable  effect  on  the  electron  density  time  dependence. 

The  measured  twilight  electron  densities  at  60  km  around  midnight 
4/5  August,  shown  in  Figure  11.4,  were  analyzed  using  a  lumped  parameter 
approach  to  find  the  time  dependence  of  detachment.  Such  an  analysis  was 
necessary  to  remove  the  effect  of  the  changing  ion  pair  production  rate 
during  this  time  period.  The  results  in  Figure  11.9  do  tend  to  show  a 
slow  rise  of  detachment  at  sunrise  which  occurs  at  approximately  the  same 
times  as  the  0  atom  rise  in  Figure  11.8.  The  computed  detachment  shows  a 
much  smaller  hysteresis  effect  than  do  the  0  atoms  in  the  1-D  calculation. 

11-19 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

lOCNHltD  tflltlltt  A  JMCI  COMPANY  INC 


LMSC/D767238 


lN3IDIdd30D  1N3WHDV13Q 


11-20 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 


•  ~  -■*-  -  ■* - -A 


£ 


Figure  11.9  Detachment  coefficients  at  60  km  altitude  calculated  from  the 
data  of  Figure  11. A.  The  sunrise  curve  is  also  shown  mirrored 
around  midnight. 
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These  results  at  60  km  based  on  measurements  and  modeling  of  the 
August  1972  event  suggest  that  the  twilight  behavior  of  electron  density 
may  be  dominated  by  the  minor  species,  0  atoms  and  (^(^A)  molecules,  rather 
than  by  photodetachment  which  is  prominent  or  dominant  in  current  negative 
ion  models.  If  this  concept  can  be  verified  by  further  study  of  all  the 
available  experimental  evidence,  reconsideration  of  present  ion  chemistry 
models  will  be  necessary. 
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12.0  THE  MASS  DISTRIBUTION  OF  IONS  IN  THE  LOWER  STK4TOSPHERE  AND  UPPER 
TROPOSPHERE 


The  mass  distribution  of  ions  is  needed  to  determine  the  mobilities  and 
hence  the  conductivities  in  the  stratosphere.  During  undisturbed  conditions, 
the  ionic  mobilities  and  conductivities  can  and  have  been  measured  directly. 
Unfortunately,  as  pointed  out  by  Meyerott,  Reagan  and  Joiner  (1980)  some 
of  the  mobilities  and  conductivities  measured  above  10  km  on  rocket  platforms 
maybe  in  error  due  to  the  high  electric  fields  employed  in  the  measurements. 
The  masses  derived  from  the  mobilities  measured  on  rocket  platforms  appear 
too  small.  The  small  masses  may  be  the  result  of  fragmentation  of  larger 
ions  broken  up  during  the  sampling  process.  Hence,  for  quiet  conditions 
improved  values  of  the  ion  masses  or  mobilities  are  required  to  provide 
an  adequate  measure  of  the  conductivities  in  the  stratosphere. 

During  disturbed  conditions,  the  mass  distribution  may  change  depending 
on  the  intensity  of  the  disturbance.  A  higher  ion  production  rate  leads 
to  higher  ion  densities.  Since  the  recombination  rate  varies  as  the 
square  of  the  ion  density,  the  ions  have  shorter  lifetimes  during  intense 
disturbances  anc  may  not  have  time  to  build  the  larger  complex  ions  that 
may  be  present  in  quiet  times.  A  measurement  of  the  mass  distribution 
can  form  a  basis  for  the  theoretical  modeling  of  the  ion  mass  distributions 
that  may  permit  one  to  predict  the  ionic  mass  distributions  as  a  function 
of  the  ion  production  rates. 

The  ionic  masses  have  been  measured  in  the  atmosphere  at  altitudes 
above  35  km  with  the  use  of  mass  spectrometers.  Ionic  masses  as  high  as 
■v  300  amu,  near  the  upper  limit  of  the  spectrometers,  have  been  detected. 
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No  insitu  mass  spectrometer  measurements  of  ionic  masses  have  been  made 
below  'v*  35  km  altitude.  Mass  distributions  can  be  inferred  at  lower  alti¬ 
tudes  from  measured  mobility  distributions  and  a  knowledge  of  the  mobility- 
mass  dispersion  relationship.  In  this  section  the  mass  distributions 
determined  from  the  measurement  of  mobility  distributions  made  in  the 
laboratory  by  Hoppel  (1970)  and  at  ^  5  km  by  Hoppel  and  Kraakevik  (1965) 
will  be  compared  with  the  mass  spectrometer  measurements  made  at  *v<  35  km 
discussed  by  Meyerott,  Reagan  and  Joiner  (1980).  Evidence  will  be 
presented  for  the  existance  of  an  appreciable  fraction  of  the  ionic  mass 
distribution  with  masses  greater  than  those  measured  with  the  mass 
spectrometers . 

If  g (m)  dm  is  the  fractional  number  of  ions  having  masses  between  m 
and  n  +  dm  and  f (k)  dk  is  the  corresponding  number  of  ions  of  having  mobil¬ 
ities  between  k  and  k  +  dk,  then 

g(m)  dm  =  f (k)  dk  (1) 


Consequently, 

g(m)  =  f(k)  (dk/dm).  (3) 

The  quantity  dk/dm  can  be  derived  from  the  mass-mobility  relationship 
given  by  Kilpatric  (1971),  reproduced  in  Figure  2  of  Meyerott,  Reagan  and 
Joiner  (1980) 

The  mobility  distribution  of  atmospheric  ions  in  a  laboratory 
simulation  has  been  measured  by  Hoppel  (1970).  The  mobility  distribution 
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of  ions  in  the  atmosphere  above  the  exchange  layer  at  altitudes  of  3.05, 
4.27  and  5.18  km  has  been  measured  by  Hoppel  and  Kraakevik  (1965). 

In  Figure  12.1  are  shown  the  mass  distributions  derived  from  the  mobility 
distributions  of  positive  and  negative  ions  measured  by  Hoppel  in  the 
laboratory  with  dry  air  and  the  negative  ions  in  moist  air.  It  can 
be  seen  in  Figure  12.1  that  the  mass  distribution  of  the  negative  ions  is 
shifted  to  the  higher  masses  as  the  water  content  of  the  air  is  increased. 

In  Figure  12.2  is  shown  the  fraction  of  the  ions  having  masses  greater  than 

m, 

G(m)  =  /  g(m)  dm  (4) 

J  m 

plotted  against  m.  It  can  be  seen  in  Figure  12.2  that  in  dry  laboratory  air 
1/2  of  the  negative  ions  have  masses  greater  than  70  amu  while  those 
in  moist  air  have  masses  greater  than  130  amu.  The  mass  distribution  of 
the  positive  ions  in  dry  air  is  weighted  toward  higher  masses  than  the 
negative  ions  with  1/2  of  the  ions  having  masses  greater  than  190  AMU. 
According  to  Hoppel,  the  addition  of  moisture  to  the  dry  air  produced  no 
appreciable  change  in  the  positive  ion  distribution. 

Laboratory  studies  of  the  nature  of  atmospheric  ions  have  been  reported 
by  Mohnen,  Del  Santo,  and  Kadlecek  (1980),  In  their  experiments  ions 
produced  in  an  ion  source  are  conducted  through  a  drift  tube  at  atmospheric 
pressure  and  the  masses  and  mobilities  of  the  emerging  ions  are  measured. 

An  electric  quadrupole  mass  spectrometer  is  used  to  identify  the  ions 
whose  mobilities  are  measured.  In  hydrocarbon  free  air  with  a  ^0  content 
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of  *  2  ppm,  the  predominant  negative  ion  is  0^  '  mass  68  amu. 

CO^,  60  amu,  C0^  ‘  (^0) ,  94  amu,  C0^  ’  (H20)^,  78  amu,  C0^,  76  amu,  0 ^  '  (H^O) , 
86  amu,  C0^  *  (820)2,  96  amu,  and  0 ^  50  amu,  in  order  of  decreasing 

concentrations  were  also  present.  With  the  addition  of  a  SO2  concentration 
of  0.14  ppm,  the  predominant  negative  ion  becomes  SO^  ,  96  amu,  with 
SO,.,  112  amu,  and  SO^  '  (^0) ,  114  amu,  also  present. 

The  laboratory  air  used  by  Hoppel  probably  contained  SO2  NO,  0^  and 
perhaps  other  impurities  in  concentration  of  parts  per  million.  Hence, 
the  negative  ion  mass  distribution  might  be  expected  to  be  broader  than 
that  reported  by  Mohnen  et  al.  (1980).  However,  it  is  interesting  to  note 
that  the  negative  ion  distribution  for  low  humidity  air  reported  by  Hoppel 
peaks  between  40  and  60  AMU  which  is  consistent  with  that  measured  by 
Mohnen  et  al.  for  pure  air,  60  to  68  amu.  This  agreement  is  reasonably 
good  considering  the  relative  crudeness  of  the  mobility  mass  determination 
and  the  fact  that  the  mobility-mass  dispersion  curve  used  may  be  less 
accurate  for  masses  below  'v  100  amu  as  can  be  seen  from  the  scatter  in 
the  points  in  Figure  2  of  Meyerott  et  al.  (1980). 

The  positive  ion  mass  distribution  derived  from  the  mobility 
distribution  measured  by  Hoppel  and  Kraakevik  (1965)  in  the  atmosphere  at 
5.18  km  is  also  shown  in  Figure  12.1.  It  can  be  seen  that  the  mass  distribution 
is  weighted  more  toward  the  higher  masses  than  those  of  the  laboratory 
mass  distributions.  From  Figure  12.2  it  can  be  seen  that  1/2  of  the  ions 
have  masses  greater  than  440  amu  and  that  20%  of  the  ions  have  masses 
greater  than  800  amu.  The  mobility  distribution  of  the  negative  ions 
differs  from  that  of  the  positive  ions  by  'v  10%  and  hence  is  not  shown. 
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The  mass  distributions  derived  from  the  mobility  distributions  of 
Hoppel  and  Kraakevik  (1963)  at  5.18  km  show  that  a  large  fraction  of  the 
ions  have  masses  in  excess  of  those  measured  with  mass  spectrometers  at 
'V  35  km  altitude  by  Arnold  et  al.  (1977),  Arijs  et  al.  (1978)  and  Arnold 

and  Henschen  (1978).  While  the  average  mobility  measured  by  Hoppel  and 

2  -1  -1 

Kraakevik  at  5.15  km  corresponds  to  a  reduced  mobility  of  1.14  cm  V  s 

which  would  correspond  to  an  ionic  mass  of  300  amu,  the  average  mass 

derived  from  the  mass  distribution  function  in  Figure  1  is  ^  600  amu. 

While  it  must  be  pointed  out  that  the  mass  distribution  function  for  the 

high  masses  derived  from  the  portion  of  the  mobility  distribution  function 

which  is  poorly  determined  and  where  dk/dm  is  large,  mobility  data  still 

suggests  the  presence  of  an  appreciable  fraction  of  the  positive  ions  at 

5.18  km  with  masses  in  excess  of  300  amu. 

The  presence  of  an  appreciable  fraction  of  the  positive  ions  with 

masses  in  excess  of  300  amu  at  5.18  km  is  not  inconsistent  with  the  mass 

spectrometer  measurements  made  at  ^  35  km  since  the  mass  spectrometers  were 

limited  to  masses  <_  328  amu.  One  might  also  expect  the  mass  distribution 

to  be  weighted  more  toward  the  higher  masses  at  5  km  than  at  35  km  due 

to  the  higher  density  and,  perhaps,  higher  moisture  and  impurity  content 

of  the  air.  However,  the  positive  ion  mobility  measurements  of  Rieckert 

2  -1 

(1971)  at  20  and  30  km  indicated  a  mean  reduced  mobility  of  *v  0.5  cm  V  s 
which  would  correspond  to  a  mass  of  •v  1000  amu.  One  is  forced  to  conclude 
that  there  is  a  strong  possibility  that  a  considerable  fraction  of  the  ions 
up  to  35  km  have  masses  higher  than  measured  by  the  mass  spectrometers. 
Simultaneous  measurements  on  the  same  (balloon)  platform  of  the  mobility 
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distribution  with  Gerdien  condensers  and  the  mass  distribution  with  mass 
spectrometers  would  be  desirable. 
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APPENDIX  A 
SPORADIC  E-LAYERS 

The  term  "sporadic  E-layer"  (or  Eg)  is  used  to  describe  a  phenomenon  of 

increases  in  the  ion  and  electron  densities  occurring  within  narrow  height 

regions  (a  few  km  thick  or  less)  at  E  region  heights  in  the  earth's  atmosphere 

(typically  between  100  and  130  km).  Figure  A.l  shows  the  altitude  profile  ior 

an  E„  at  temperate  latitude  (Smith  and  Miller,  1980).  In  this  case  the  maxi- 
s 

mum  densities  increased  by  nearly  two  orders  of  magnitude  over  the  ambient 
background  in  a  thin  2-3  km  slab. 

There  are  different  classes  of  Eg  probably  having  different  origins;  for 
the  middle  latitude  Eg  the  preferred  explanation  assumes  concentration  of 
metallic  ions  in  one  or  more  narrow  layers  by  wind  shears.  At  high  latitudes 
energetic  particle  precipitation  may  be  the  source  (Jones  and  Saksena, 

1980).  Whitehead  (1970)  reviews  the  extensive  Eg  literature  and  identifies 
three  main  types,  equatorial,  temperate,  and  auroral  zone  Eg.  Figures  A. 2  - 
A. 5  (Smith,  1957)  show  some  of  the  relevant  characteristics  for  the  latter  two  types 
It  will  be  noted  that  in  the  auroral  zones  Es  is  typically  an  evening  or  night¬ 
time  phenomenon  since  the  occurrence  frequency  drops  markedly  during  the  hours 
around  midday.  At  temperate  latitudes,  however,  maximum  occurrence  is  ob¬ 
served  around  the  summer  solstice  with  predominance  of  daytime  occurrences 
(Figure  A. 3). 

However,  from  the  point  of  investigating  possible  effects  of  Eg  on  ELF 
propagation  important  information  seems  to  be  lacking  or  is  insufficient  for 
realistic  calculations,  e.g.  little  is  known  of  the  spatial  extent  of  the  Eg 
layers,  some  may  be  1000  km  or  more  in  extent  (DeGregorio  et  al.,  1962), 
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whereas  others  may  be  100  km  or  less  in  horizontal  extent  (Whitehead,  1970). 
Quantitative  information  on  homogeniety  of  ionization,  vertical  and  horizontal 
movements  of  the  Eg,  geographical  positions  at  given  times,  time  duration, 
etc.  is  furthermore  limited. 
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ELECTRON  DENSITY  (cm'5) 

Figure  A.l  Double-peaked  sporadic  E-layer  observed  near  midnight 
(Smith  and  Miller,  1980). 
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Figure  A. 2  Mean  annual  Eg  occurrence  in  the  auroral  zone  fEc  >  5  MHz: 
1948-1954 
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APPENDIX  B 

THE  EFFECTS  OF  ENERGETIC-PARTICLE  PRECIPITATION  ON  THE  GLOBAL  ATMOSPHERIC 
ELECTRICAL  CIRCUIT 

Solar  energetic  particle  precipitation  during  large  solar  particle  events 
(SPE)  greatly  enhances  the  conductivity  of  the  atmosphere  over  extended  spa¬ 
tial  regions.  In  addition,  for  the  larger  events  the  electrical  current 
carried  by  the  solar  protons  and  electrons  in  the  polar  caps  greatly  exceeds 
the  normal  air-earth  current  that  flows  between  the  ionosphere  and  the  ground 
at  high  altitudes.  In  order  to  model  the  effects  of  the  SPE's  on  the  earth's 
electrical  system,  we  note  that  the  precipitation  is  reasonably  uniform  over 
the  polar  caps  whose  dimensions  are  several  thousand  km.  Since  the  important 
perturbation  occurs  over  altitudes  below  100  km.  the  problem  is  essentially  a 
plane  problem.  We  have  used  the  SPE's  of  August  1972  as  an  example  of  large 
SPE's  since  the  particle  precipitation  that  occurred  during  that  period  is 
better  documented  than  that  for  most  events. 

The  electric  field  in  the  fair  weather  regions  (no  thunder  storms  present) 
can  be  calculated  by  assuming  that  the  electrical  current  density  is  indepen¬ 
dent  of  altitude.  This  assumes  quasi-steady  state  conditions  and  that  cur¬ 
rents  are  only  a  function  of  altitude.  The  air-earth  current  density,  JA£,  is 
the  sum  of  the  current  density  carried  by  the  precipitating  protons,  Jp,  the 
current  density  carried  by  the  precipitating  electrons,  Jg,  and  the  ohmic 
current  density,  Jq 
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We  have  used  the  usual  convention  that  a  downward  flowing  positive  charge  is  a 
positive  current.  Since  the  current  carried  by  the  precipitating  particles  is 
zero  near  the  earth,  J^g  is  equal  to  the  air-earth  current  near  the  earth. 

From  Equation  (1)  is  given  by 


J 


52 


AE 


+  J 

e 


(2) 


The  electric  field  E(H)  at  any  altitude,  H,  is  given  in  terras  of  the  conduc¬ 
tivity,  o(H)  by 


E(H)  =  V£(H)  (3) 

o(H)  depends  on  the  particle  deposition  rate  and  the  mobility  of  the  ions. 

The  conductivity  in  the  polar  cap  changes  during  a  SPE  due  to  the  greatly 
increased  deposition  of  high  energy  protons  from  the  SPE  and  also  due  to  the 
reduction  of  galactic  cosmic  rays  known  as  the  Forbush  decrease.  Changes  in 
the  flux  of  galactic  protons  at  the  ground  are  usually  measured  by  neutron 
monitors.  Figure  B.l  is  a  neutron  monitor  record  taken  at  the  South  Pole 
during  the  period  3-11  August  1972  from  Poraerantz  and  Duggal  (1973).  It  can 
be  seen  in  Figure  B.l  that  the  first  Forbush  decrease  (FD)  started  at 
~0300  UT  on  August  4.  At  ~1200  UT  on  August  4  a  sharp  increase  in  the  neutron 
counting  rate  occurred  and  continued  until  ~2000  UT.  This  so-called  Ground 
Level  Event  (GLE)  occurred  at  the  same  time  as  the  10-hour  intense  solar 
proton  precipitation.  Following  this  GLE,  a  second  FD  occurred.  This  second 
FD  lasted  for  ~4  days.  At  ~1500  UT  August  7,  a  second  GLE  occurred.  An 
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increase  in  solar  proton  precipitation  was  also  recorded  by  the  satellites  at 
the  time  of  this  second  GLE. 

The  ion  production  rates  due  to  the  precipitation  of  solar  cosmic  rays 
and  galactic  cosmic  rays  are  usually  determined  by  different  techniques  due  to 
the  difference  in  the  spectrum  of  the  precipitating  protons.  The  solar  pro¬ 
tons  have  energies  ranging  from  ~1  MeV  to  ~1  BeV  while  the  galactic  primary 
protons  have  energies  from  ~500  MeV  to  ~100  BeV.  The  principal  energy  loss 
mechanism  for  protons  with  energies  less  than  ~500  MeV  is  ionization  loss, 
while  for  protons  with  energies  greater  than  ~500  MeV  the  principal  energy 
loss  mechanism  is  nuclear  collisons.  Protons  with  energies  of  ~1  BeV  or  less 
produce  low  energy  protons  and  neutrons  whose  ranges  are  small  compared  to  the 
primary  protons.  Hence,  the  energy  lost  by  the  primaries  in  nuclear  colli¬ 
sions  can  be  considered  to  produce  ionization  at  the  point  of  the  nuclear 
collision.  At  proton  energies  greater  than  1  BeV,  mesons  are  produced  as  well 
as  high  energy  protons  and  neutrons.  The  mesons  have  ranges  much  larger  than 
the  ranges  of  the  primary  protons  so  that  the  energy  loss  by  the  primary 
protons  cannot  be  considered  as  a  local  deposition.  In  fact,  only  about  1%  of 
the  ionization  at  sea  level  is  produced  by  protons.  Mosts  of  the  ionization 
is  produced  by  mesons. 

Ion  production  rates  due  to  galactic  cosmic  rays  have  been  extensively 
measured  as  a  function  of  altitude,  latitude,  and  solar  cycle  (Neher,  1971  and 
references  therein).  The  method  used  to  determine  the  ion  production  rates 
during  a  Forbush  decrease  is  to  reduce  the  empirical  ion  production  rates 
determined  for  a  given  altitude,  latitude,  and  time  of  the  solar  cycle  by  the 
square  root  of  the  fractional  reduction  in  counting  rate  of  the  neutrons 
monitors  from  the  pre-SPE  value. 
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The  ion  production  rates  for  the  solar  protons  have  been  determined  in 
the  past  at  Lockheed  PARL  from  the  proton  spectrum  determined  from  satellite 
measurements  augmented  at  the  high  energy  end  by  balloon  measurements  (Reagan 
et  al.,  1980)  through  the  use  of  the  PROTON  deposition  code.  This  spectrum 
which  has  been  used  for  the  previous  studies  at  LMSC  extends  from  ~1  MeV  to 
~400  MeV.  For  the  evaluation  of  the  effects  of  peak  proton  precipitation 
during  the  period  ~1200  UT  to  ~2000  UT  August  3,  it  appeared  necessary  to 
extend  the  solar  proton  spectrum  to  higher  energies  so  that  there  would  be  a 
better  overlap  with  the  ion  production  rtes  produced  by  the  galactic  protons 
at  altitudes  below  20  km.  The  evidence  for  significant  proton  flux  at  ener¬ 
gies  greater  than  400  MeV  is  evidenced  by  the  first  GLE  shown  in  the  neutron 
monitor  data  in  Figure  B.l. 

The  solar  proton  spectrum  for  energies  greater  than  ~400  MeV  during  the 
first  GLE  on  August  4,  1972  has  been  determined  by  Lockhwood  et  al.  (1975) 
from  analysis  of  the  word-wide  neutron  monitor  data  taken  at  that  time.  We 
have  used  the  spectrum  determined  by  Lockhwood  to  extend  our  previously  deter¬ 
mined  solar  proton  spectrum  for  energies  greater  than  400  MeV.  The  combined 
solar  proton  spectrum  for  energies  greater  than  100  MeV  for  1508  UT  August  4 
is  shown  in  Figure  B.2.  Also  shown  in  Figure  B. 2  for  comparison,  is  the 
galactic  cosmic  ray  proton  spectrum  for  1972  given  by  Fulks  (1975).  We  have 
used  this  solar  proton  spectrum  as  representative  of  the  solar  proton  spectrum 
for  the  entire  ~10  hour  peak  irradiation  period. 

The  extension  of  the  solar  proton  spectrum  to  include  protons  with  ener¬ 
gies  above  400  MeV  required  some  modification  of  the  PROTON  code  to  take 
nuclear  collisions  into  account.  This  was  done  by  multiplying  the  proton  flux 
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Figure  B.2  Solar  Proton  spectra  for  energies  greater  than  100  MeV  for 
period  ^  1200  to  'v  2000  UT  A  August  1972.  Galactic  cosmic 
ray  proton  spectrum  for  1972. 
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of  energy  E  making  an  angle  0  with  the  vertical  at  a  height  h,  F(E9),  by  exp 

(-/j“  p(h)dh/Acos9),  where  p(h)  ■  atmospheric  density  at  height  h  and  A  =  mean 

2 

free  path  for  nuclear  reactions.  The  value  of  A  =  205  gm/cm  as  suggested  by 
Peters  (1958)  for  protons  of  energy  ~400  MeV  was  used.  The  energy  lost  per 

3 

cm  through  nuclear  collisions  was  assumed  to  be  deposited  locally  and  was 
added  to  that  lost  through  ionization  at  a  rate  of  35  eV  per  ion  pair. 

The  PROTON  deposition  code  was  also  modified  so  that  it  would  calculate 
the  proton  current  across  a  horizontal  unit  area  at  any  altitude  as  well  as 
the  energy  deposited  at  that  altitude.  In  Figure  B  .3  is  a  plot  of  the  proton 
current  density  divided  by  3.70  x  10  amp.m  as  a  function  of  altitude  at 
1508  UT  August  4,  1972,  the  most  intense  period  of  that  SPE.  The  nominal  air- 
earth  current  density  in  the  polar  cap  is  3.70  x  10  amp/m  .  It  can  be  seen 
that  the  proton  current  density  becomes  comparable  to  the  air-earth  current  at 
~30  km  and  is  ~50  times  larger  than  the  nominal  air-earth  current  at  ~100 
km.  Also  shown  in  Figure  8.3  is  the  electron  current  density  that  was  estimated 
from  the  measurements  made  at  Chatanika,  Alaska.  While  the  proton  current 
density  is  fairly  uniform  over  the  polar  cap,  the  electron  current  exhibits 
considerably  spatial  variation  due  to  auroral  activity.  Since  Chatanika, 

Alaska  is  in  the  auroral  zone,  the  electron  current  density  shown  in  Figure 
B.3  is  probably  an  upper  limit  to  the  current  density  in  the  polar  cap. 

The  ion  conductivity,  can  be  calculated  from  the  positive  and  nega¬ 

tive  ion  concentrations,  N+  and  N_,  and  the  respective  mobilities,  K+  and  K_, 
by 

ot  -  eN+K+  +  eN_K_  (4) 


B-7 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

lOCKMlIO  MISSIIIS  *  IMCI  COMPANY.  INC 


CAP  PRECIPITATION  CURRENTS 


LMSC/D767238 


CM 


I 

O 


X 

IS. 


o 


2 


B-8 


LOCKHEEO  PALO  ALTO  RESEARCH  LABORATORY 


l  O  C  «  N  (  I  0 


Mismn 


S  f  k  c  * 


c  o 


I  N  C 


Figure  B.3  Polar  cap  currents  carried  by  precipitating  solar  protons  and  electrons 
expressed  as  the  ratio  of  a  typical  air-earth  current. 
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In  the  lower  atmosphere,  below  50  km,  N+  2  N_,  and  K+  ~  K_  so  that 

oi  -  2eNK  (5) 

When  the  Ion  production  rates,  Q,  are  known,  N  is  given  in  terms  of  the  pro¬ 
duction  rates  and  the  ion-ion  recombination  coefficient,  by 

N  «*  (Q/a,  ,)1/2  (6) 

1  9  1 

Hence 

oL  -  2eK(Q/ai  1)1/2  (7) 

In  Figure  B.4  is  shown  the  calculated  ion  conductivity  in  the  polar  cap  for 
quiet  periods  for  1972.  The  production  rates  were  assumed  to  be  due  to  cosmic 
rays  for  that  period  of  the  solar  cycle.  Conductivities  for  three  different 
values  of  the  mobilities  corresponding  to  ions  masses  of  100,  200,  and  1000 
AMU  are  shown.  Also  shown  are  conductivity  measurements  by  Kraakevik  (1961) 
for  altitudes  from  0  to  6  km  and  by  Tanaka  et  al.  (1977)  from  20  to  35  km. 
the  derived  conductivities  of  Tanaka  et  al.  were  obtained  from  measurements  of 
the  air-earth  current  and  electric  fields  in  the  polar  caps,  with  o  = 

The  conductivity  measurements  of  Kraakevik  were  made  over  Greenland  in  1954. 
This  is  approximately  the  same  relative  period  in  the  solar  cycle  as  the 
August  1972  SPE.  The  measurements  by  Tanaka  et  al.  were  made  over  the  South 
Pole  in  November  -  December  1972. 
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It  can  be  seen  in  Figure  B.4  that  the  conductivities  derived  from  the 
measurements  of  Kraakevik  fit  the  calculated  conductivity  curve  for  an  ion 
mass  of  100  AMU  while  the  measurements  of  Tanaka  et  al.  fall  between  the 
curves  for  a  mass  of  200  and  1000  AMU.  As  discussed  by  Meyerott,  Reagan,  and 
Joiner  (1980)  the  low  value  of  the  conductivity  at  ~30  km  may  be  due  either  to 
ion  densities  at  this  altitude  which  are  lower  than  that  predicted  by  the 
theory  or  because  the  masses  of  the  ions  at  this  altitude  are  greater  than  200 
AMU. 

The  quiet  time  ionic  conductivity,  o,  ,  for  August  1972  was  taken  to  be 
the  curve  shown  in  Figure  B.4  which  fits  the  conductivities  derived  from  the 
measurements  of  Kraakevik  from  0  to  6  km  and  those  derived  from  the  measure¬ 
ments  ot  Tanaka  et  al.  at  ~30  kra.  The  ionic  conductivity  disturbed  period 
1508  UT  August  4,  a ^  was  obtained  from  the  undisturbed  ionic  conductivity 
by 


l,d 


=  a.  (Q./Q  ) 

i,u  d  u 


1/2 


(8) 


This  approach  assumes  that  the  mobilities  and  the  recombination  coefficients 
during  the  disturbance  are  the  same  as  those  during  the  undisturbed  times. 

The  disturbed  ion  production  rates,  Q^,  were  taken  to  be  the  sum  of  the  pro¬ 
duction  rates  by  galactic  cosmic  rays,  qR,  and  the  production  rates,  Qgpg, 
calculated  from  the  measured  proton  spectrum  from  the  modified  PROTON  code  as 
described  previously. 


"  Qd,CR+ 


SPE 
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As  can  be  seen  from  the  neutron  monitor  data  shown  in  Figure  1  the  magnitude 
of  the  Forbush  decrease  during  the  peak,  of  the  event  was  8%;  hence,  is 
given  from  the  undisturbed  cosmic  ray  production  rates,  by 


Q 


d,CR 


0.92 


hi.CR 


(10) 


In  Figure  B.5  is  shown  the  conductivity  in  the  polar  cap  for  altitudes  up  to 
35  km  for  the  disturbed  period  1508  UT  A  August  1972  and  for  the  undisturbed 
periods  of  1972.  Figure  B.6  shows  the  conductivities  extended  from  35  km  to 
90  km.  It  can  be  seen  that  the  disturbed  conductivity  is  slightly  smaller 
than  the  undisturbed  below  ~12km  and  is  much  larger  than  the  disturbed  for 
altitudes  above  20  km. 

The  electric  field  for  the  disturbed  and  undisturbed  periods  can  be 
calculated  from  Equation  (3).  For  the  disturbed  period.  Equation  (3)  becomes 


E(H) 


AE  i,d 


-  J  /o  ,  +  J  la.  . 
p  i,d  e  i,d 


(ID 


In  Figure  B.7  the  vertical  electric  field  for  the  disturbed  and  undisturbed 
case  is  shown  as  a  function  of  altitude.  It  can  be  seen  that  the  disturbed 
electric  field  becomes  negative  from  ~30  km  to  ~55  km.  This  is  due  to  the 
fact  that,  as  is  shown  in  Figure  B.  3  f  the  proton  current  is  larger  than  the 
air-earth  current  for  altitudes  above  ~30  km.  At  ~55  km  the  electron  pre¬ 
cipitation  current  exceeds  the  proton  current  and  hence  the  electric  field 
becomes  positive  above  that  altitude.  Electric  field  measurements  have  been 
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Figure  B.6  Calculated  undisturbed  and  SPE  disturbed  conductivities  on  4  August  1972  at 

Chatanika,  Alaska  at  0508  local  time  (1508  UT)  during  tbe  main  peak  of  the  event. 
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made  at  30  kta  by  Holzworth  and  Mozer  (1979)  during  the  5  August  1972  SPE. 

They  measured  a  very  small  electric  field  at  30  km.  Considering  the  error  In 
their  measurements  It  Is  not  possible  to  tell  whether  the  field  was  positive 
or  negative. 

The  Implications  of  the  change  in  conductivities  and  the  large  precipi¬ 
tation  currents  in  the  polar  cap  on  the  potential  of  the  ionosphere  can  be 
inferred  by  calculating  the  ionospheric  potential  inside  the  polar  cap  and 
just  outside  the  polar  cap  assuming  the  same  air-earth  current  at  both  loca¬ 
tions.  Inside  the  polar  cap  the  electric  field  is  the  disturbed  field  shown 
in  Figure  B.7.  Outside  the  polar  cap  the  conductivities  are  the  undisturbed 
conductivities  reduced  by  the  Forbush  decrease  of  ~0.4%.  In  Figure  B.8  is 
shown  the  difference  in  potential  at  1508  UT  as  a  function  of  altitude  between 
positions  inside  (V^)  and  outside  (VQ)  the  polar  caps,  for  an  air-earth  cur- 
rent  density  of  3.70  x  10  arap/m  .  Since  the  conductivity  both  inside  and 
outside  the  polar  cap  is  changed  by  the  same  amount  up  to  ~10  km  due  to  the 
Forbush  decrease,  the  potential  up  to  this  altitude  remains  the  same  at  the 
two  positions.  At  30  km,  however,  the  atmospheric  potential  at  the  two  loca¬ 
tions  differs  by  almost  10  kilovolts  due  to  the  presence  of  the  solar  particle 
currents  and  the  enhanced  conductivities  in  the  polar  caps  but  not  outside  the 
polar  caps. 

A  paper  entitled  "The  Effects  of  Energetic-Particle  Precipitation  on  the 
Global  Atmospheric  Electric  Circuit"  based  on  the  work  in  this  report  was 
presented  at  the  VI6*1  International  Conference  on  Atmospheric  Electricity  held 
in  Manchester,  England,  28  July  to  1  August,  1980. 
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Figure  B.8  Difference  in  the  earth's  potential  at  1508  UT,  4  August 

1972  as  a  function  of  altitude  between  positions  inside  (V^) 
and  outside  (VQ)  the  polar  caps. 
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